SUM-INTEGRAL INTERPOLATORS AND THE EULER-MACLAURIN FORMULA 

FOR POLYTOPES 



STAVROS GAROUFALIDIS AND JAMES POMMERSHEIM 

Abstract. A local lattice point counting formula, and more generally a local Euler-Maclaurin formula 
follow by comparing two natural families of meromorphic functions on the dual of a rational vector space V, 
namely the family of exponential sums (S) and the family of exponential integrals (I) parametrized by the set 
of rational polytopes in V. The paper introduces the notion of an interpolator between these two families 
of meromorphic functions. We prove that every rigid complement map in V gives rise to an effectively 
computable Si-interpolator (and a local Euler-MacLaurin formula), an IS-interpolator (and a reverse local 
Euler-MacLaurin formula) and an IS°-interpolator. Rigid complement maps can be constructed by choosing 
an inner product on V or by choosing a complete flag in V. The corresponding interpolators generalize the 
work of Berline-Vergne, Pommersheim-Thomas, and Morelli. 
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1. Introduction 



1.1. What is a local lattice-point counting formula? The relationship between sums and integrals 
has been of interest to mathematicians since the ancient Greeks. The classical Euler-Maclaurin formula, 
discovered in the first half of the eighteenth century shortly after the development of modern calculus, may 
be viewed as a relationship between the sum of a function over the lattice points in a one-dimensional 
polytope with integer vertices and the integral of the function over the polytope. One naturally asks the 
same question in higher dimensions: Given a polytope P in an n-dimensional space V equipped with an 
n-dimensional lattice A and a function / on V, can one express the sum of / over the lattice points in P in 
terms of the integral / over P? In such a formula, one would expect a main term involving the integral of 
/ over P as well as correction terms involving the integrals of / over the proper faces of F C P. 

In the simplest case, let us suppose that / is a constant function. Then, the question becomes that of 
expressing the number of lattice points in P in terms of volume of P and the volumes of the faces F of P. 
If P is 2-dimensional, the celebrated Pick 's formula [Pi] 



expresses the number #(P) of lattice points in a convex lattice polygon in terms of its area A = vol(P), and 
the number of lattice points b = vol(dP) of its boundary. For example, we have: 



Here and throughout, we follow the usual convention that all volumes vol(P) of faces are normalized so that 
a lattice basis of AflL, where L is the linear space parallel to F, has volume 1. Unfortunately, Pick's formula 
is not local, due to the presence of the term 1. Here, locality means that for each face F of P, the formula 
contains a term that is the volume vol(P) multiplied by a coefficient that depends only on the supporting 
cone Supp(P, F) to P along F. The supporting cone is defined as the union of rays whose endpoint is in 
F and which remain in P for a positive distance. See Section 12 . 1 1 below for a precise definition. McMullen 
[McMj proved the existence of local lattice point counting formulas. More precisely, he proved the existence 
of a function /i from rational, convex cones to rational numbers such that for any integral polytope P, the 
number #(P) of lattice points in P is given by 



Here the sum is taken over all faces F of P, and vol(P) denotes the volume of the face F. 

Part of the difficulty constructing and computing a function \i that satisfies ([T]) is that /i is far from unique. 
The second author and Thomas gave an explicit construction of a rational valued function /i satisfying ([1} , 
given a fixed complement map, a notion introduced in Thomas's thesis (cf. |Thj ) . A complement map 
is a systematic choice of complements of linear subspaces of a vector space; see Section 12.41 for a precise 
definition. Also note that all of the complement maps in this paper will be rigid. Two natural ways to get a 
rigid complement map are to choose: (1) an inner product, or (2) a complete flag. Given an inner product 
or a complete flag, |PT) construct effectively a map fi that satisfies Equation ([I}. For the triangle depicted 
above, the Pommersheim- Thomas values of fi at the supporting cones to the vertices for an inner products 
or complete flags are as follows: 



#(p) = A + h + i 





(1) 




F 



_3_ 
10 




d% — Adi d2-\-1d\ 
6d 2 (d 2 -2di) 




1 
4 



9 d\ +3dic£2+d<2 



lldj — 7d\d2-\-d\ 
12di(2di-da) 



20 12did 2 



SUM-INTEGRAL INTERPOLATORS AND THE EULER-MACLAURIN FORMULA FOR POLYTOPES 



3 



The values on the left arise from the standard inner product on Z 2 , according to |PT[ Cor.l], and may be 
computed by multiplying out the Todd polynomial in the ring presentation given in |PT[ Prop. 2]. Readers 
wishing more details on this should consult Section 15.21 The values in the figure on the right are those 
arising from the choice of the complete flag in V* = Z 2 whose 1-dimensional subspace is spanned by the 
point (dijcfo)- These values can be computed using the method outlined in |Mo| p. 198]. They can also be 
computed by multiplying out the Todd polynomials using |PT| Thm.3]. One feature of this construction is 
that if one chooses a complement map arising from a complete flag, one recovers exactly the lattice point 
formulas of Morelli |Moj , who gave a function /i that satisfies Equation ((T|) and takes values on the field of 
rational functions on a Grassmannian. For example, in the triangle on the right above, the values shown 
are rational functions on the Grassmannian of 1-dimensional subspaces of V*. The construction of |PT] is 
based on the theory of toric varieties, and gives an answer to a question of Danilov about the existence of a 
local expression for the Todd class of a toric variety. 

1.2. What is a local Euler-MacLaurin summation formula? Returning to the Euler-Maclaurin ques- 
tion, Berline and Vergne constructed in |BVj an explicit local Euler-MacLaurin formula for the sum of a 
polynomial function / over the lattice points Pn A of an n-dimensional polytope P in a rational vector space 
V with lattice A. The Berline- Vergne formula has the form 

(2) £ m = J2 1 D ^ F )-f 

igPnA F Jf 

where the sum is over the set of faces F of P, and D(P, F) an infinite-order constant-coefficient differential 
operator D(P,F) that depends only on the supporting cone Supp(P, F). Equation ([2]) is a generalization of 
(fT]). Indeed, if D(P,F) satisfies <j2j) , and we define ^(Supp(P, F)) to be the constant term of D(P,F), then 
the local lattice point formula ([T]) holds. 

As in McMullen's case, the infinite order differential operators D(P,F) are not uniquely determined by 
([21). The construction of Berline- Vergne requires an inner product on the vector space V, and their results 
apply to rational polytopes, rather than just integral polytopes. Their construction associates to each cone 
K in V a meromorphic function [i(K) on the dual space that is regular at the origin. For cones of dimension 
at most 2, the value of this function at (the constant term of the operator D(P,F)) coincides with the fi 
constructed in jPT] for the special case of complement maps arising from an inner product. In this way, for 
the example triangle above, the Berline- Vergne construction recovers the the numbers j,-§q,jq arising out 
of the inner product case of the Pommersheim-Thomas construction. This fact that this coincidence holds 
is new to this paper; see Theorem [TJ properties (5) and (6). 

1.3. An informal presentation of the results of this paper. The local Euler-MacLaurin formula @ is 
a consequence of a relationship between the integral and the sum of an exponential function over a polytope. 
In this paper, we introduce the concept of an interpolator between the families of exponential sums (S) and 
exponential integrals (J) over rational polytopes in a rational vector space V. We now discuss our results 
about interpolators informally, leaving the precise definitions and statements for Sections and [31 

If P is a rational polytope, or more generally a rational polyhedron (region, not necessarily compact, 
defined by linear inequalities with rational coefficients) in a rational vector space V, one can associate to 
P two important meromorphic functions, the exponential sum S(P) G M(V*) and the exponential integral 
I(P) £ M.(V*) where M.(V*) is the algebra of meromorphic functions on the dual space V* ®C The values 
of these functions at a point £ £ V* ® C are given by: 

(3) S(P)(0= Yl ^ I(P)(Z)= f e^dx 

xePnA Jp 

provided \e^' x ^ | is summable (resp. integrable) over P. Here, the integral is taken with respect to the relative 
Lebesgue measure on aff(P), normalized so that a basis of A n W, where W is the linear subspace parallel 
to the affine span of P, has volume 1. The fact that Equations define meromorphic functions, as well as 
the precise characterization and properties of the functions S and /, is essentially the content of Lawrence's 
theorem [La], reviewed in Section l2~2l below. 
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The local Euler-Macluarin formula follows in a straighforward manner from a formula of the following 
shape for a rational polyhedron P 



where [i is a function on the set of cones in V with values in Ai(V*). We call such a function an SI- 
intcrpolator, or simply an interpolator. The main result of this paper (TheoremQJ states that a complement 
map on the vector space V gives rise in a natural way to an effectively computable Si-interpolator on V, 
and hence a local Euler-Maclaurin formula of the form [5] see Theorem [5J 

Thus, in particular, an inner product on V or a complete flag in V gives rise to a local Euler-Maclaurin 
formula. These interpolators have interesting connections with previous results. For complement maps 
arising from the choice of an inner product, one recovers the /i constructed in |B Vj . In addition, we show 
that the values of these functions fi at are given by the functions \i in [PTj for cones of dimension at most 
2, and conjecture that these coincidences hold in all dimensions. In the case of complete flags, one obtains 
interpolators /x that are entirely new, as can be seen by calculations in dimension 2. In this case, the values 
of the constant term fj,(K)(0) coincide with values constructed by Morclli for cones of dimension at most 2, 
and conjecturally in all dimensions. One can also vary the chosen flag, thus associating to each cone K a 
meromorphic function fi(K) on F1(V*) x V* , where F1(V*) is the complete flag variety of V. This function 
n{K) is naturally defined and effectively computable independent of any choices; see Theorem[3l In Morclli's 
work, he constructs what amounts to the constant term //(i"T)(0) of this meromorphic function. Working 
dimension-by-dimension, Morelli chooses to view this term as a rational function on a Grassmannian, rather 
than on the entire flag variety. 

In Section 13. 2\ we show that a complement map on V also leads naturally to the construction of a local 
IS-interpolator, expressing the exponential integral I(P) in terms of the exponential sums S(F), for F a 
face of P; see Theorem 2J This construction, which works only in the case of integral polyhedra, allows one 
to obtain a reverse local Euler-MacLaurin formula expressing the integral of a polynomial function over a 
polytope in terms of the sums of the function over lattice points in the various faces of the polytope; see 
Theorem [6l 

Finally, we show how a complement map yields an IS°-interpolator, where we use S° to denote the sum 
over interior lattice points of a polytope. Part of the interest in IS°-interpolators lies in an observation of 
Morelli, who asked essentially (in the language of the present paper) if the same function can simultaneously 
serve as the constant term of an Si-interpolator and as the constant term of an IS°-interpolator on the dual 
space. He observed that the constant terms that he constructs do exactly this in dimensions at most 4. We 
give a proof of these coincidences in dimensions 1 and 2, and extend these results to the case of complement 
maps arising from inner products. Our arguments are essentially geometric, in contrast to Morelli's rather 
algebraic argument. (See Sections 13.31 and l6l for a complete discussion.) 

Section [5] contains explicit computations of the functions (J,(K) for some cones K of dimensions at most 
2. We compute these functions for complement maps arising both from inner products and from complete 
flags. In Section 15721 we exhibit Equation @] in the case of two triangles including the triangle depicted above, 
for both the inner product and complete flag cases, matching the constant terms with those constructed by 
Pommershcim-Thomas and Morelli shown in the above figure. 



2.1. Polytopes. We will use standard terminology for rational polytopes, polyhedra, cones and faces, fol- 
lowing for example |Fu] , We will fix a finite dimensional lattice A, isomorphic to Z fc for some k £ N, and 
consider the rational vector space V = A ®z Q- All polyhedra will be rational and convex, and all cones will 
be rational, polyhedral and may be affine. A polytope is a compact polyhedron. If P is a polyhedron and F 
is a face of P, then the tangent cone Tan(P, F) is defined by picking a point x in the relative interior of F 
and looking at all the directions that one can go and stay in P: 



(4) 




F 



2. Polytopes, exponential sums/integrals and interpolators 



Tan(P, F) = {v e V\x + ev e P for small e > 0} 
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This cone is independent of the choice of x and contains the origin. Let aff(P) denote the affine span of F, 
i.e., the smallest affine subspace of V that contains F. Let lin(P) denote the linear subspace of V parallel to 
aff(P). lin(P) is the maximal linear subspace of V contained in Tan(P, F). The supporting cone Supp(P, F) 
is the tangent cone translated back to its original position: 

Supp(P, F) = Tan(P, F) + x 

for x in the relative interior of F. 




2.2. Exponential sums and integrals. Let V* = Hohiq(V, Q) denote the dual vector space to V, and 
(•, •) : V* x V — > Q denote the natural evalutation pairing on V . By a meromorphic function on V* we 
will mean a meromorphic function on the complexified dual space V£, where Vc — V <8>q C denotes the 
complcxification of V. Let Ai(V*) denote the algebra of meromorphic functions on V* and A4 r (V*) denote 
the subalgebra of meromorphic functions regular (i.e., analytic) at the origin. Note that if W is a rational 
subspace of V, then A' = A n W is a full rank lattice in W which naturally induces a relative Lebesgue 
measure on W, normalized so that a basis of A' has volume 1. Likewise, every rational affine subspace of V 
has a natural measure induced by the lattice. 

Given a polyhedron P in a rational vector space V with lattice A as above, there are two important 
meromorphic functions S(P) and I(P) with rational coefficients in V* , namely the generating function of 
exponential sums and exponential integrals. 

These functions are uniquely characterized by the following properties. 

(Al) If P contains a straight line, then S(P) = I{P) = 0. 

(A2) S (resp. I) is a valuation (resp. a solid valuation). That is, if the characteristic functions of a 
family of polyhedra satisfy a relation J2i r iX(Pi) — 0> then the functions S(Pi) satisfy the relation 
J2i TiS(Pi) = (resp. restrict the sum to those Pi that do not lie in a proper affine subspace of V.) 

(A3) For every s €V and every £ £ V* we have 

(5) I(s + P)(0 = e<^I(P)(0, 
and 

(6) S( S + Pm = e^S(P)({;), seA. 

(A4) If £ S V* is such that |e^^| is integrable (resp. absolutely summable) over P, then 

(7) I(P)(0 - / e^dm P (x), S(P)(0 = eiM 

•' p xePnA 

where dmp denotes the relative Lebesgue measure on aff(P). 
The existence of the functions I and S is by no means obvious, and is essentially the content of Lawrence's 
theorem; [La] . (Al)-(A4) determine the value of I and S on a simplicial cone as follows. 

(A5) If K = Cone(vi, . . . , Vk) is a simplicial cone generated by k independent vectors v±, . . . , Vk, (where 

k = dim(y) in the case of S) and □(i>i, . . . , Ufe) denote the parallelepiped generated by v\, . . . ,Vk, 

we have: 

<8) HKm = (-i)* yoli ^ f -f\ smo-( E Hnr^y 
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Since any polyhedron can be subdivided (virtually or not) into a union of simplicial cones, together with 
Equation © above, one obtains an algorithm to compute I(P) and S(P). The complexity of such an 
algorithm is discussed in [Bvj . The next theorem which computes I(P) and S(P) of a polyhedron in 
terms of the tangent cones of its vertices was obtained by Brion (using toric varieties) and Lawrence (using 
combinatorics) independently. 

Proposition 2.1. [Brl ILaj Let P be a polyhedron in V. Then. 

(9) I(P)= /(Supp(P,t;)), S(P)= S(Supp(P,v)) 

VEV(P) 'UGV(P) 

where the sums are over the set V(P) of vertices of P. 

2.3. Si-Interpolators. So far, our discussion of the meromorphic functions S(P) and I(P) was parallel but 
independent: our formulas did not mix S and I. It is natural to ask for an interpolation between S(P) and 
I(P). Indeed, if P is a polytope, then letting £ = reduces to a classic problem: expressing the number of 
lattice points in P in terms of the volumes of the faces of P. Let C(V) denote the set of cones of V. 

Definition 2.2. An Si-interpolator (or simply, interpolator) on V is map 

fj, : C{V) — ■* M{V*) 

such that for any rational polyhedron P in V we have: 

(10) S(P)= M(Supp(P,F))J(F) 

where the sum is over the set F(P) of all faces of P. An interpolator is regular if fi takes values in A4 r (V*). 

For the next definition, recall that if W is a quotient of V with projection tt : V — > W, then there is a 
natural map C(W) — > C(V) given by K n> tt^ 1 ^). Moreover, W* is naturally a subspace of V* , and there 
is a restriction map M(V*) — > M(W*). 

Definition 2.3. An Si-interpolator /i on V is hereditary if for every rational quotient W of V, the composition 
ft : C(W) — > M(W*) given by the following diagram 

C{W) C(V) 



M(W*) M(V*) 

is an interpolator on W. 

The next lemma reduces the checking of the interpolator equation (|10p to the case of cones. 

Lemma 2.4. fi is an interpolator if and only if it satisfies (|10|) for all cones P in V . 

Proof. This follows from Proposition 12. II Indeed, we have: 

S(P) = Yl S(Supp(P,vj) 
vev(P) 

= E E /i(Supp(Supp(P,«),F'))/(i^). 

vev(P) F'eJ r (Su PP (p : u)) 

Since every face of Supp(Supp(P, v), F') is a face of Supp(P, F) for some face F of P, it follows that 

S(P)= Yl m(Su PP (P,P)) Y /(Supp(F,«))= E M(Supp(P,P))7(P). 
FeJ r (p) vev(F) FeF(P) 

The result follows. □ 
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2.4. Rigid Complement maps. A rigid complement map gives us a systematic way to extend a function 
on a linear subspace of a vector space to the entire vector space. This notion was introduced in the thesis 
of Thomas |Thj and used in |PTj . where it is shown that a choice of complement map leads naturally to 
a local formula for the number of lattice points in a polytope and a local formula for the Todd class of a 
toric variety. In our paper, we will only use the notion of a rigid complement map. The curious reader may 
consult [FTl HE] for the defmit ion of the general complement maps. 

Definition 2.5. A rigid complement map on V* is a map from a collection C* of linear subspaces in V* 
to the set of linear subspaces of V* satisfying two properties: 

(a) For every U G C*, $>(U) is complementary to U; i.e., *(J7) n U = {0} and + U = V*. 

(b) If Ut C U 2 C V* with L7 lf f/ 2 e C*, then *(C/ 2 ) C *(Z7i). 

There are two easy ways to construct rigid complement maps: (a) an inner product Q on V*, and (b) a 
complete flag L on V". 

Lemma 2.6. (a) An inner product Q on V* defines a complement map with domain the set of all 
subspaces of V* such that 

where t/^ denotes the subspace of V* perpendicular to U C V* under the inner product Q. 

(b) A complete flag L = (La, . . . , L n ) on V* satisfying Lq = {0} C Li C - • • C L ra = V* defines a complement 

map on V* by 

* (^0 = £n-dim(t/) 

for any linear subspace U of V* which is complementary to L n -dim([/) ■ This complement map is only defined 
generically, i.e., for subspaces that meet the flag generically. 

Our main theorem, Theorem [1] below, requires cones that are generic with respect to a rigid complement 
map tf>. Let us define what generic means. Fix a rigid complement map on V*. 

Definition 2.7. (a) A quotient W of V is ^-generic if W* € C*, where W* C V*. 

(b) A subspace U of V is ^-generic if V/U is ^-generic. 

(c) A cone K of V is ^-generic if lin(F) is ^-generic for every face F of K. 

(d) A polyhedron P of V is ^-generic if for every face F of P, Supp(P, F) is ^-generic. 

The next lemma states that the notion of a complement map on V* is hereditary, i.e., a complement map 
on V* gives rise to a complement map on every W* , where is a ^-generic rational quotient of V. 

Lemma 2.8. If ^ is a rigid complement map on V* and W is a ^-generic rational quotient of V , then 
there is a natural complement map $ on W* defined by 

for all U C W* such that ^(U) CI W* is complementary to U . 

The next lemma gives the promised extension of functions on a linear subspace W* to ones on the entire 
space V*. 

Lemma 2.9. If ^ is a rigid complement map on V* and W is a \& -generic quotient ofV, then the decom- 
position 

V* — W* © *(W*) 

defines a unique linear projection map 

7T : V* — > W* 

that annihilates ^(W*) and is the identity on W* . This allows any function on W* to be extended to V* . 

We now come to a main definition. Fix a rigid complement map ^ on V* . Let C (V) denote the set of 
^-generic cones of V. 

Definition 2.10. A ^-compatible interpolator /i on V is a function [i : '(V) — > Ai(V*) satisfying 
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(a) Equation (|10p for all ^-generic polytopes P, 

(b) /i is ^-hereditary, i.e., for every ^-generic quotient W of V, the composition p, : C(W) 
given by the following diagram 



M(W*) 



(11) 



C*(W) 



M{W*) 



C*(V) 



M{V* 



satisfies Equation (|T0|) for all ^-generic polytopes P in W. 
(c) Moreover, the following diagram commutes: 



C*(W) 



C*(V) 



(12) 



A4(T/T) — A<(V*) 



where n is given in Lemma 12.9 



3. Statement of the results 

3.1. SI- interpolators and local Euler-MacLaurin formula. Our main result is that a rigid complement map 
^ on V* determines algorithmically a unique ^-compatible interpolator /i* on V. This generalizes results 
of |BV1 IMol IPT| . The proof uses the techniques of |BVj . which were a motivation and inspiration for us. 

Theorem 1. (a) If V is a rational vector space and \1/ is a rigid complement map on V* , there is a unique 
^/-compatible interpolator /i*. 

(b) In addition, /i* satisfies the following properties on the collection of vE'-generic cones on V: 

(1) (Additivity) If the characteristic functions of a finite collection of cones Ki with vertex v 6 V satisfy 
the relation J^. r,x(iQ) = 0, then the functions satisfy the relation J^. {Ki) = . 

(2) (Lattice Invariance) If v G A, then /x*(u = (^(K). 

(3) (Isometry Equivariance) If g is a lattice-preserving linear automorphism of V and g* is the inverse 
transpose, then ^ (g(K))(g*(Z)) = M*(*0(0- 

(4) (Regularity) ^i* is regular at £ = 0. 

Moreover, ^ generalizes previous results of [Moj and |BV] : 

(5) (Constant Term) Let K be a top dimensional cone in V with vertex at zero, and suppose that V 
has dimension at most 2. Then for arising from an inner product or a complete flag, the constant 
term of ^ (K) agrees with the Todd class coefficient /i(AT v ) of the dual cone K v that depends on 
ty, constructed in |PT| Cor.l]. In particular, for "J/ coming from a complete flag, the constant term 
of (K) gives Morelli's formula |Mo) . 

(6) (Inner Product) For ^> coming from an inner product Q* on V* , the function /j, (K) agrees with 
n(K) constructed in [BVj . where K is the image of K in V/V(K), where V(K) denotes the linear 
subspace parallel to the largest affine subspace contained in K . 

There are two corollaries of Theorem [TJ a local version of the Euler-MacLaurin formula for polytopes 
(promised in the introduction), and a meromorphic function on Fl(y*) x V* associated to a pointed cone K 
in V, where F1(V™) is the variety of complete flags of V. 

To formulate these results, fix a rigid complement map ^ on V , a ^P-generic rational convex polyhedron 
P and a face F of P. Consider the Taylor series expansion of /x*(Supp(P, F)){£,) and the corresponding 
differential operator (of infinite order, with constant coefficients): 



(13) 



D*(P,F)= M *(Supp(P,F))(d x 
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Theorem 2. Let \& be a complement map on V* , let P be a \& -generic rational polytope in V , and let h(x) 
be a polynomial function on V . With the above notation, we have 

(14) h &= E / D*(P,F)-h. 

xePnA FeT(P) 

In particular, if h(x) — 1 then 

(15) #(P) = J2 M*(Supp(P,P))(0)volOF) 

where #(-P) (resp. vol(P)) denotes the number of lattice points (resp. the volume) of P. 

For concrete examples, see Section [3] Equation (I15[) computes the number of lattice points of a polytope 
as a weighted sum of the volume of its faces. 

Suppose now that \1/ comes from a flag L. Letting L vary, we obtain the following. 

Theorem 3. Given a pointed cone K in V there exists an effectively computable meromorphic function 
l-i(K) in F\(V*) x V* which is regular in a Zariski open subset of F1(V*) x {0}. 

For an algorithmic computation of ^(K) for cones K of dimension at most 2, see Section [5] 
The next conjecture identifies the constant term of a ^-compatible interpolator with the one of |PT| . We 
post it as a conjecture for now, and hope to discuss it in a future publication. 

Conjecture 1. Part (5) of Theorem^ holds in all dimensions. 

Question 1. Is there any relation between the rational functions n(K) in Fl(y*) x V* and the quantum 
cohomology ring on ¥\(V*)1 

3.2. IS-interpolators and reverse Euler-MacLaurin formula. The definition of the SI- interpolator leads in a 
natural way to the notion of an IS-interpolator, which may be a useful notion to numerical approximations 
of integrals by sums. The next notion requires us to restrict attention to lattice polyhedra, i.e., polyhedra 
whose vertices are points of the lattice A. Below, a lattice cone K in V is a rational cone such that Uf)A ^ 0, 
where U is the largest affine subspace contained in K. 

Definition 3.1. An IS-interpolator is map A, from the set of all lattice cones in V to Ai(V*) such that for 
any lattice polyhedron P in V we have: 

(16) I(P)= A(Su PP (P,F))S(F) 

where the sum is over all faces F of P. An IS-interpolator is regular if A takes values in A4 r (V*). 

Hereditary and ^-compatible IS-interpolators are defined in an analogous way to Definition 12.31 and 12. 101 
The statement and proof of Theorem Q] holds with minor modification for IS-interpolators. 

Theorem 4. (a) If V is a rational vector space and \1/ is a rigid complement map on V* , there is a unique 
^/-compatible IS-interpolator A . 

(b) In addition, satisfies properties (2)-(5) of Theorem [1] and the following version of additivity: 

(1) (Additivity) If K is a cone that does not contain a linear subspace and is subdivided into a finite 
union of cones K i: then A* (AT) = E l: dim(K 1 )=di m (K) 

One can obtain directly a formula for A* in terms of the interpolator /x* of Theorem [T] as follows, by 
observing that /i*(Supp(if, K)) = 1. 

Theorem 5. If K is a positive (resp. zero) dimensional lattice cone then 

(17) Yl A*(F)^(Su PP (X,F)) = 0, (resp.) A*(tf) = 1. 
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The interpolator A* gives a reverse Euler-MacLaurin summation formula. Fix a rigid complement map 
W on V, a ^-generic lattice polyhedron P and a face F of P. Consider the Taylor series expansion of 
A*(Supp(P, P))(£) and the corresponding differential operator (of infinite order, with constant coefficients): 

(18) A*(P ) F) = A*(Supp(P,F))(a a: ) 

Theorem 6. With the above assumptions, for every lattice polytope P and every polynomial function h{x) 
on V we have 

(19) / h(x)= E (&*(P,F)-h)(x). 



FeJ r (p) xeFnA 



In particular, if h(x) — 1 then 



(20) vol(P) = A*(Su PP (P,P))(0)#(P) 

where #(P) denotes the number of lattice points of F. 



Equation (|20|) computes the volume of a lattice polytope as a weighted sum of the number of lattice points 
of its faces. A formula of this type was first written down by Morelli, |Mo[ Eqn.(5)]. 



3.3. IS°-interpolators and Morelli's work. When ^ comes from a complete flag, Equations (fT"5)) and 
(|20|) are similar with some results [Mo( Eqn.(5),Eqn.(6)] of Morelli. To explain this, let us introduce the 
variant S° of the exponential sum function defined by: 

(21) S°{P) = S{P°) 

where P° denotes the interior of a polyhedron P. P° is not a polyhedron itself, however it is a virtual sum 
of polyhedra. Then, we can talk about IS°-interpolators v on lattice cones. Theorems [4] and [5] have the 
following analogue. 

Theorem 7. (a) If V is a rational vector space and \& is a rigid complement map on V* , there is a unique 
^-compatible IS°-interpolator i/*. In addition, 

(b) (Additivity) If K is a cone that does not contain a linear subspace and is subdivided into a finite union 
of cones K it then v* [K) = E I:d im(K I )=dim(K) 

(c) satisfies properties (2)-(5) of Theorem[T] 

(d) If if is a lattice cone then 

(22) J2 v*(F)»*(Supp(K,F)) = l. 

Fix a lattice polytope P and a face F of P. If 4" comes from a complete flag L in V* as in Lemma [2?6l and 
we vary the flag, we obtain two rational functions i^(Supp(P, F j) and /x(Supp(P, F j) on F1(V*) X V* which are 
regular at a Zariski open subset of F1(V*) X 0. Thus, we can consider the constant terms ^*(Supp(P, P))(0) 
and /i*(Supp(P, P))(0) which are rational functions on F1(V*). 

On the other hand, Morelli constructs a pair of rational functions Vk^pP) and piki'&FP) (where k is 
the dimension of lin(Supp(P, F))) that appear respectively in Equations (5) and (6) of |Mo) (unfortunately, 
Morelli denotes the two functions with the same notation, although he clearly recognizes that they are 
distinct functions). Morelli's construction uses K-theory, localization and the Bott residue theorem. 

The following proposition is a re-expression of the flag case of Theorem [TJ part (5). 

Proposition 3.2. With the above notations, for polytopes P of dimension at most 2, we have: 

(23) i/(Supp(P,P))(0) = u k {M), MSu P p(P,F))(0) = n k {p F P). 
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In [Mol p. 191] and also in [Mo| Thm.5], Morelli observes that 



(24) »{K)(Q) = K^ V )(0) 

for lattice cones K of dimension 4 or less, and that Equation (l24l) fails in dimensions more than 4. This 
strange coincidence is easy to explain from our point of view, for cones of dimension at most 2. The rational 
functions v(K) and /i(AT v ) are distinct, in fact they are functions on distinct spaces (F1(V*) x V* and 
F1(V) x V, respectively). Their constant terms coincidentally agree in small dimensions. For a discussion, 
see Section 



4. Proofs 

4.1. Proof of Theorem [TJ Fix a rigid complement map f on a rational vector space V with lattice A, and 
consider a ^-compatible interpolator /i* and a Vt-generic rational polytope P. We will show by induction 
on the dimension of P the existence and uniqueness of fi (P). All cones in this section will be assumed 
to be \l/-generic. By Lemma [2.41 it suffices to assume that P = K is a cone in V. To avoid confusion, we 
denote the map /j, : C (V) — > A4(V*) by fj,y and for every ^-generic quotient W of V, we denote the map 
P : C*{W) — ► M{W*) by (m w . 

If K is pointed with vertex v, then we can single out the contribution from the O-dimensional face of K 
in Equation (|10|) , and together with hereditary property of Equation ([TT]) , it follows that 



(25) M *(if)(0 = e-<W £ e«'*>- £ 4 /lin{F) (Supp(K,F))(n(0) [ e^dm F (x) 

\xGKnA F,dim(F)>0 F 

where Supp(AT, F) denotes the image of Supp(P, F) in the quotient space V/lin(F), and it : V* — > 
(Vyiin(F))* is the projection using the decomposition V* = (V/]m.(F))* ® *((V/lin(F))*) given by the 
rigid complement map v & 

On the other hand, if if is a non-pointed cone in V, let V(K) denote the linear subspace of V parallel to 
the largest affine subspace contained in K, and consider the image K of K in V/V(K), which is a pointed 
cone of dimension strictly less than the dimension of K . Equations (jlip and ()12ll imply that 

(26) $(K)(Q=4 /V(K) (K)(0 

where £ e V*, it : V* — > (V/V(K))* is the projection using the decomposition V* = (V/V{K))* © 
^>((V/V(K))*) given by the rigid complement map ^ and ^ € (V/V(K))* is the image of £ under the above 
projection. 

Equations (|25f and ([26| uniquely define fiy from for dim(V4 7 ) < dim(l/). On the other hand, when 
V is a O-dimensional space, we have /^*({0}) = 1. This uniquely determines [iy. An explicit computation 
of /^(K) for cones K of dimension 0, 1 and 2 is given in Section [SJ 

Additivity of /z* follows from the above inductive definition of /i* and the additivity (property (A2)) of 
the exponential sum and exponential integrals. This is presented in detail in the proof of |BV[ Prop. 15]. 

Remark 4.1. In fact, Berline-Vergne in |BV| work with the transverse cone T(i^, F) of a face F of K defined 
to be the image of Supp(A', F) under the projection map V — > V/lm(F). T(K, F) is always a pointed cone. 
In fact this gives a 1-1 correspondence 



(27) {cones in V} O {pointed cones in quotients of V} 

Using this correspondence, and an inner product Q on V, Berline-Vergne inductively construct \fi defined 
on the set of pointed cones on the quotients of V . For visual reasons, we prefer to work with the supporting 
cones Supp(AT, F) rather than the transverse cones T(K,F). 
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Going back to the proof of Theorem [TJ the lattice invariance and the isometry equivariance of /n* follows 
easily by induction and the invariance of the / and S functions. 

The regularity of the meromorphic function /x* (K) at zero follows by induction and a residue calculation, 
discussed in detail in |BV[ Prop. 18], using the following fact about the exponential sum and integral functions. 

(A6) If if is a cone in V with primitive integral generators V\, . . . ,v^, then 

k k 

(28) llfoviWK)®, Hit^IiKm 

i=l z=l 

are regular functions and the residues of the meromorphic functions S(K) and I(K) along the 
hyperplane v± — satisfy the equations 

(29) Res Vl (S(K)) = -S(K), Res Vl {I(K)) = -I(K). 

(A6) follows in turn by additivity and Equation ((5J) which evaluates the / and S functions on simplicial 
cones. 

We defer the proof of Property (5) of Theorem [1] until the example section (Section [5]). See the proofs 
immediately following Examples 15.61 and 15.71 

Finally property (6) of Theorem Q] follows directly from |BV[ En.4] . 

4.2. Proof of Theorem [2j Fix a polynomial function h(x) on V and a polytope P on V. Consider the 
constant coefficient differential operator h(d^) which acts on the function £ i-^ e^'^ by: 

(30) h(d 6 )e^' x) = h{x)e^' x) . 
The definition (J?} of S(P) and the above gives that 

h(d £ )S(P) = Hx)e {i ' x) 
xePnA 

Let ev denote the evaluation at £ = 0. It follows that 

(31) (evoh(d i ))S(P)= £ h(x). 

xePnA 

On the other hand, Theorem [1] gives that: 

S(P)(0 = £ M *(T(P,F))(0/(F)(0 

= E / l**(T(P,F))(Oe {tx) dmp(x) 

F£j r (P) 

E / D*{P,F)e^ x Um P {x) 
Fer(P) F 

where the last equality follows from the definition ([13")) of the differential operator D*(P, F). Applying the 
differential operator h(d^) to both sides it follows that 

h(d ( )I(P)(0= E I D^{P,F)h{x)e^ x) dm P {x). 

FET(P) Jf 

Evaluating at £ = 0, it follows that 
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(32) (evoh(ds))S(P)= £ f D*(P,F)-h. 

Equations (|3 1 1) and (|32|) complete the proof of Theorem [2] 

4.3. Proof of Theorem |4j The proof of Theorem [1] applies verbatim, with the following observation, which 
explains the need for lattice polyhedra. 

If K is a rational pointed cone with vertex v for which (| 16[) applies, then one of the terms in (|16[) is 
F = {v}. In that case, T(A, {v}) = K and the corresponding term in (p~6|) is v(K)S({v}). Now, (A3) and 
(A5) imply that 




if v e A 
otherwise. 



Thus, when K is not a lattice cone, then we cannot solve for v(K). In the case of Si-interpolators, as 
discussed in Section [2~4l the corresponding term of (TTOt was /i(K)I({v}) and (A3) and (A5) imply that 

/({«})(£) = e«'">. 

4.4. Proof of Theorems \S\ and [3 Theorem [5] follows from the fact that the composition of an SI- 
intcrpolator with an IS-interpolator is an SS-interpolator, and the fact that an SS-interpolator is unique. 
More precisely, fix an Si-interpolator /i and an IS-interpolator A, and consider a lattice cone K. Then, on 
the one hand we have 

S{K)= J2 M(Supp(A',F'))/(F'). 

f 'eF(K) 

On the other hand for every face F' of K, we have 

I(F')= J2 HSnpp(F',F))S(F) 

FeF(F') 

Substituting the second equation into the first, it follows that 

S(K)= J2 ^(Sup P (A,F'))A(Supp(F',^))5(F). 

F,F': FCF' 

This motivates the following definition. Consider the function 

Supp(F, F) ^ (A o M )(Supp(P, F)) := ^ A*(Supp(A, F'))A(Supp(F', F)). 

F': FCF' 

It is easy to see that Xo/j, is an S'S'-interpolator. Moreover, if /x and A are ^/-compatible, so is their composition. 
On the other hand, there is a unique VP-compatible S'S'-interpolator. Thus (A o ^)(Supp(P, F)) = Oresp. 1 
for F 7^ P (resp. F = P). This concludes the proof of Theorem [5] 

Theorem [7] follows by an analogous computation, using the fact that the map (y o /i) is a ^-compatible 
SS°-intepolator, the uniqueness of such interpolators, and the fact that the map that sends every cone to 1 
is a ^/-compatible SS°-intepolator. 

5. Computations and Examples 

5.1. Computation of fi for low dimensional cones. As was mentioned earlier, the meromorphic func- 
tions /z*(Af) are effectively computable given a complement map In this section we illustrate this, by 
explicitly computing /i*(A) for cones K of dimension at most 2. 

Proposition 5.1. Suppose that K = {0} in the vector space V — {0}. Then fx(K) is the constant function 
1, independent of ^ . 

Proof. S{K) = I{K) = 1 forces fi(K) = 1. □ 
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If {vi, . . . , Vk} is a collection of vectors in V, let (y\, . . . , Vf.) = Q + «i + ■ ■ • + denote the cone spanned 

by those vectors, where Q + is the set of non-negative rational numbers. 

Proposition 5.2. Now suppose that V is one- dimensional and K = Cone(u) is a ray in V generated by a 
primitive vector v G A. Then, independent of ^ , we have 

pi{K){0 = B((tv)), 

where 

(33) B{z) = Th + \-\-T2 + ^-^ + --- 

is the generating series of the Bernoulli numbers. 

Proof. Use the property of interpolators 

S{K) = M (Supp(K, K))I(K) + ^(Supp(K, 0))/(0), 

/i(Supp(AT, K)) — 1 by the previous proposition and compatibility under quotients. In addition, we have 
S{K){£) = 1/(1 - e«) and I{K){£) = The result follows. □ 

For the following proposition, recall that if K is a subset of V, then we define the dual K y by 

(34) K v = {w € V* | (w, v) > 0, for all v £ K}. 

Also, if A is a lattice in V, then A* = {w G V* \ (w, v) G Z, for all v G A} is the dual lattice in V*. Below, 
we will denote by fi L (resp. fjP) the ^/-compatible Si-interpolator where comes from a complete flag L on 
V (resp. an inner product QonV), as in Lemma 12.61 

Proposition 5.3. Suppose V is two-dimensional, and K <zV is a half-plane with boundary a line U through 
the origin. 





Let p G A* be the primitive generator of the ray K v . Suppose the 'J is defined by a complete flag L in V* , 
and let c be a generator of LI, where L\ is the one- dimensional subspace of the flag L (the only nontrivial 
subspace in the flag.) Then 

Now suppose W is defined by an inner product Q on V . Then 



H L (K)(0 = B 



p Q (K)(0 = B 



Q{Lp) 
Q(p,p) 



Proof. We use the ^-compatibility together with our one-dimensional formula above. To do so, we must 
compute the projection from ir : V* — > W* where W = V/U . This projection is computed using the 
decomposition V* — W* "i>*(W*) given by the complement map 1 J r . Suppose the 1-dimensional subspace 
^(W*) C V* is generated by d. Then for £ G V* , find scalars uji,oj 2 such that 

(35) £ = Ldip + UJ2d, 

Then from ^-compatibility and the one-dimensional formula, we obtain 
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For ^ coming from a complete flag, we pair both sides of Equation (|35[) with the generator c. Since (d, c) = 0, 
we find 

Wi - 7 r, 

and the formula for p. L (K) follows. For \I> coming from the inner product Q, we note that Q(d, p) = 0. Thus 
Equation [35] yields 

U1 = Q&£) 

Q{p,py 

which completes the proof. □ 

Using the defining property of interpolators, we immediately get the following proposition, an explicit 
expression for p, of a nonsingular 2-dimensional cone. 

Proposition 5.4. Suppose V is two-dimensional and K = Cone(ui, V2) where v±, v 2 form a basis of A C V. 
Let Fi — Cone(fi) and F 2 = Cone(i>2), and let p\-,p% G A* denote the primitive normals to F± and F%. Then 
for \& coming from a complete flag L in V* , we have 

r ,<n u l (K] = 1 1 1 1 r{ fc c > V 1 n( fc c > 

( ) Mlj (l-e<^>)(l-e<^>) (Z,v 1 )(t,v 2 ) + (Z, Vl ) {(p 1 ,c)) + (Z,v 2 ) {(p 2 ,c) 

For coming from an inner product Q on V* , we have 

M :: Q (K] 1 1 I 1 p( Q & Pl) \ I 1 E f Q ^ P2) 

[ ' ^ [ ' (l-e«.M)(i- e «.«»>) (Z,vi)&v 2 ) + <£, Wl ) \Q(pi,Pi)J (t,v 2 ) \Q(p 2 ,p 2 ) 

Remark 5.5. Replacing 1/(1 — e x ) with B(x) — 1/x, it follows that Equations (j3"6"]l and ([5T|) can be written 
in the form: 



(f,W2) V \Q{P2,P2)J 

5.2. Examples. In this section we explicitly compute some low dimensional examples and match them with 
the ones given by Morelli, Pommersheim-Thomas and Berline-Vergne. 

Example 5.6. Consider the triangle P with vertices wo = (0,0), V\ — (1,0), v 2 = (0,1) in the lattice A = 
1? C V = Q 2 , where V* is also identified with Z 2 . 




Let ^ be the complement map coming from the complete flag L = (Lq, Li,L 2 ) in V* such that L\ = 
Span{(di, e^)}- To compute p L of the supporting cones, we use lattice invariance to translate them to the 
origin. In other words, consider the tangent cones Ki = Tan(P, V4) for i = 0, 1, 2, to P at V4. Specifically we 
have 

K a = Cone((l, 0), (0, 1)), K x = Cone((-l, 1), (-1, 0)), K 2 = Cone((0, -1), (1, -1)). 



Then for £ = (£1,6) S V* , we have 



16 



STAVROS GAROUFALIDIS AND JAMES POMMERSHEIM 



(i L (Ko)(£) = T , FT ~ TT + T B (^ ~ TT^ + T B (& ~ ^) 

(1 - e«i)(l - e^) 66 6 V rf 2 J 6 V d i J 

M (KlKO = T JZTTTm =7F + TT7 7T ~ & - 3-6 + 7 T B 



(l-e&-fr)(l-e-fr) 6(6-6) 6 V di V 6-6 V di-da 
{K2m -(i-e^)d-e-&) + 6(6^6) ~ 6 B l a " *M + 6~ 6 S l^T^ 

These functions are analytic at £ = (0,0). For example, ^ l (-Kq)(6 has the Taylor expansion 



" (/fo)(0 = ( — m& — J + I24 J €l + \ rmaW\ l fl $ 



;2 



-3d?-t?V 2 (-U\-d\ 



£<= _L t 

12M\d 2 ) V 720^di ,v 
Denoting the constant term of each \i L (Ki) by ^(Ki), we find 



L/ zr \ 4 + 4+ 3rfirf 2 
^ (ifo) = 12^ 



Mo(^ 
Mo TO 



12di(di -tfe) 
5df - 5did 2 + d\ 
12d 2 (d 2 -d 1 ) 



in exact agreement with Morclli's example |Mol p. 198-199]. These rational functions can also be computed 
by multiplying the Todd polynomials using jPTl Theorem 3] In particular, these three rational functions 
sum to 1 in agreement with Equation (|15p. 

In addition, one checks by direct computation that 

J^ L (Supp(P,P))/(F) = S(P) = 1 + e? 1 +e« 2 . 

F 

in agreement with the definition of an interpolator. 

Indeed, this example provides a general proof of the coincidence of our construction and the flag case of 
[PTj for cones of dimension at most 2, and hence the coincidence with Morelli's formulas. 

Proof of QJ part (5) for Hags: In dimension 1, the [PTj values are always equal to h, independent of 
complement map. Using Proposition 15.21 we see that this coincides with the constant term of /x. For the 
two-dimensional case, we first use the additivity of [PTj and of our construction (Theorem 1, Property (1)) 
to reduce to nonsingular 2-dimensional cones. By functoriality, we may assume that K is the cone Kq from 
example I5.6[ generated by the standard basis of Z 2 . In this case, we note that the values of /Uq (Kq) given in 
Example 15.61 match the values in |Mol p. 198-199]. By |PT1 Corollary 2], these values also match the values 
in |PTj . This completes the proof. 



Example 5.7. We consider the same triangle P from Examplc l5.61 but now with a complement map ^ coming 
from an inner product Q on V*, given by the matrix 

a b 
b c 
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We can use Proposition 15.41 to compute pP(Kj) for the supporting cones Ki of P. We find 

f fi ( R o)(0 = 71 cTTn W\ -TT + T B U + + T B (& + 7& 

(l-e«)(l-e") 66 6 V a J 6 V c 

M W (^i)(0 = 7i iwTTTl — 777 + T77 7T ~ 6 + 76 + 7 7" S ' 



(l_ e &-£i)(l_e-Ci) 6 V c V 6 -6 V a + 26 + c 

q^w^ 1 1 1 , b \ 1 R /-a6-&6-&6-c6 

^(^2X6 = 7; "77~77rT7i — 777 + T77 7T ~ 6 + 76 I + 7 



(1 - eft-^)(l - e-C») 6(6-6) 6 V a 7 6~6 V a + 2b 
This time, we find the constant terms are given by 



M?(^i) 



3ac — ab — be 
12ac 

afe + 4ac + 106c + 26 2 + 5c 2 

12(ac + 26c + c 2 ) 
5a 2 + 2fc 2 + lOab + Aac + be 



12(a 2 + 2a6 + ac) 

These agree with the values of p constructed in |PT[ Cor.l] for coming from an inner product. In 
particular, the three rational functions above sum to 1. For the standard inner product on V* = Z 2 , 
(corresponding to a = c = 1, b = 0), we find 

^{K,) = \ **?(*i) = f ^{K 2 ) = \. 

For completeness and comparison, we detail how these numbers arise, in a seemingly very different way, 
out of the construction of Pommersheim-Thomas. We begin, according to |PT[ Cor.l(iv)] by considering the 
outer normal fan of P. The rays of this fan are po — (1, 1), p\ — (— 1, 0), and p 2 — (0, —1). We must then 
multiply the second-degree Todd polynomial T = J2i<j DiDj + J2 Si in the ring given in |PT[ Prop. 2]. 
In this ring, one finds D 2 = ^D Di + ^D D 2 , D\ = D Q Di, and D\ — D D 2 . Hence 

T = \d x D 2 + ^D D 2 + |d Di, 

from which one reads off the values of po(Ko) = j, po(Ki) = |, and po(K 2 ) = |. In fact, more generally, 
let K = cone(vi, v 2 ) be any nonsingular cone in V with vertex at 0, and let Q be an inner product on V* . 
Then an easy computation similar to the above can be used to compute the value of po of [PTj for this cone 
K. Indeed, if K* — cone(pi, p 2 ) is the dual cone, then using presentation of |PT| Proposition 2], one finds 
that the coefficient of D\D 2 in T is given by 



f 38 ) 1 Q(PUP2) ( 1 1 \ 1 Q*{vi,v 2 ) ( 1 1 \ 

1 ' 4 12 \Q( Pl ,pi) Q( P2 ,p 2 )J 4 12 \Q*(v ltVl ) Q*(v2,v 2 )J' 

where Q* denotes the dual inner product on V. 

Indeed, this calculation proves the inner product case of Theorem [T] part (5), as we now show. 

Proof of \J\ part (5) for inner products: As in the flag case, the [PTj values are always equal to |, 
independent of complement map, matching the constant term of the functions in Proposition 15.21 For the 
two-dimensional case, we again use additivity and functoriality to reduce to the case in which K is the cone 
Kq from example 15.61 We then note that the formula of Equation [38] for the [PTj construction matches the 
values of Pq(Kq) given in Example 15.71 This completes the proof. 

A comparison of the functions p^(Ki) of Example 15.71 with the functions pg(Ki) of Example 15.61 reveals 
that they are indeed different, confirming the fact that interpolators that come from flags are different from 
those that come from an inner product. 

Let us repeat the previous example with a different triangle. 
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Example 5.8. Consider the triangle P with vertices v = (0, 0), Vi = (2, 0), v 2 — (0, 1). 




As above, let ^ be the complement map coming from the complete flag L = (Lo, Ly, L 2 ) in V* such 
that LI = Span{(di, ^2)}- Denoting, as before, the tangent cones at the vertex Vi by Ki, we see that 
Kq = Cone((l, 0), (0, 1)), and K\ = Cone((— 1, 0), (— 2, 1)) are nonsingular cones, whose p values may be 
computed from Proposition 15.41 On the other hand, the cone K 2 — Cone((0, —1), (2, —1)) is singular. We 
can compute p L (K 2 ) by subdividing K 2 into two cones 

K 2 =C 1 UC 2) Ci = Cone((0,-l),(l,-l)), C 2 = Cone((l, -1), (1, -2)) 

using the ray p = Cone((l, —1)). The additivity of p implies that 

p L {K 2 )=p L (C 1 )+p L (C 2 )-p L (p). 



The terms on the right hand side can be computed using Propositions 15.41 and 15.21 Alternatively, one 
can compute pr{K 2 ) directly from the defining formula for interpolators. Letting F\ — Cone((0, —1)) and 
F 2 = Cone((2, —1)) denote the one dimensional faces of K 2 , we get 

p L (K 2 ) = S{K 2 ) - I{K 2 ) - M i (Supp(^ 2 ,F 1 ))/(F 1 ) - p L (Su VV (K 2 ,F 2 ))I{F 2 ) 
1 + e*-* 2 , 1 B U^ + 1 



(l-e-£2)(l-e 2 ex-e2) 6(26-6) -6 V d 2 s y 26-6 V d 2 -2d, 

The constant terms of p L (K\) are Morelli's rational functions 

L df + 3did 2 + dl 

A^o K K o) = 



12did 2 



L lldf - 7d x d 2 + d\ 



12di(2di - da) 
L d\ - 4d x d 2 + 2d\ 

M*a)- -6d 2 (2d 1 -d 2 ) 

agreeing with the numbers shown in figure in the introduction. 

Again, one can verify by direct computation that 53j=o Vai^-i) = 1 an d 

5^/i £ (Supp(P,i0)/(-F) = S(P) = 1 + eP + e 2?1 + e« 2 . 

F 

Example 5.9. We consider the same triangle P from Example l5.81 but now with a complement map coming 
from an inner product Q on V*, given by the matrix 

a b 
b c 

As in Example 15.81 the function p® of the singular cone K 2 may be computed either by subdividing, or by 
direct use of the defining property of interpolators. We find that pfi(K 2 ) is given by 

1 + , 2 , J_ B U_ b_A + ( -«6 - 266 - 66 - 2c6 " 



(l-e-« 2 )(l-e 2 ^-^) 6(26-6) -6 V a ) 26-6 V a + 46 + 4c 
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The constant terms, in agreement with the construction of [PTj . are 

o , T , s 3ac — ab — be 



\2ac 



Q ab + 5ac + 46 2 + 256c + 22c 2 



12(ac + 46c + 4c 2 ) 
Q _ 2a 2 + 8ab + lac + 2b 2 + 2bc 
Mo ( 2) " 6( a 2 + 4a6 + 46c) 
which add up to 1. For the standard inner product, these coefficients are 

Again for comparison, let us detail how these numbers arise out of the construction of Pommersheim- 
Thomas. To apply |PT[ Cor.l(iv)], we first consider the outer normal fan of P. The rays of this fan are 
Po = (1,2), pi = (—1,0), and P2 = (0,-1). A nonsingular subdivision is achieved by adding the ray 
P3 = (0, 1). As before, we multiply the Todd polynomial T — J2i<j DiDj + i Yli &l m the ring given in 
[PTl Proposition 2]. One calculates D 2 = \{2D Q D 2 - 2D D 3 ), D 2 = 0, B\ = 2D D 2 , and D 2 = -2D D 3 . 
Hence 

T = \d x D 2 + ^D Q D 2 + ^D a D 3 + jDiAs, 

Using additivity under subdivision |PT1 Cor.l(iii)], one finds pq(Kq) = j, p (Ki) — and po(K 2 ) = 
j_ , l j_ 

20 ' 4 10' 

6. Comparison with Morelli's work 

In this section we compute values of v and discuss the relations Morelli noticed between the constant 
terms of p and v. 

6.1. Computation of v for low dimensional cones. For cones of dimension 0, 1 and 2 it is easy to 
compute v (K) explicitly, just as we did for p, (K) in Section [5j The proofs in this case are similar. 

Proposition 6.1. Suppose that K — {0} in the vector space V = {0}. Then v{K) is the constant function 
1, independent of ^ . 

Proposition 6.2. Now suppose that V is one- dimensional and K — Cone(u) is a ray in V generated by a 
primitive vector v G A. Then, independent of ^ , we have 

Proposition 6.3. With the notation of Provosition [573[ we have 
and 

Proposition 6.4. With the notation of Proposition \5~^j\ we have 



<^«i)«,«a> 
and 



«.«0 V (Pi,c)J l-e«.«») V (P2,c)J (i_ e «.*i>)(l-e«>«»)) 



^ V Q(Pi,Pi)J l-e«.«»> V Q(P2,P2)J ' (l-e«-"i>)(l-e«^>) 
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6.2. Morelli's coincidence in dimensions 1 and 2. In this section we will prove Morelli's observation 
for lattice cones of dimension at most 2, using the computations of v^(K) for cones K of dimension 
at most 2 from Section 16.11 Using these propositions, we can prove the following relation between fi and v 
for cones of dimension at most 2. The following theorem, proved by Morelli in the flag case, shows that the 
coincidence of constant terms holds also in the case of inner products. 

Theorem 8. Let K be a cone with vertex at in V , and assume that the dimension of V is at most 2. Let 
^ be a complement map on V induced by an inner product on V* or a complete flag in V* and let be the 
corresponding complement map on V . Then we have 

Proof. For cones of dimension 0, the constant term of both p, and v is always equal to 1. For cones of 
dimension 1, the constant term of both /i and v is always 1/2. Now let K have dimension 2. From the 
inversion formula of Theorem it follows that 

Hence, the desired conclusion is equivalent to 

M^)(o) + M^ v )(o) = i 

Suppose that fi is induced by an inner product. The inner product gives us an identification of V with 
V* , allowing us to consider K and K v as both living in the same space V. If K = Cone(ui, V2), then 
K y = Cone(wi, w%) with m orthogonal to Wi, i = 1,2, under the inner product. Furthermore, considering 
the orthogonal cones L\ = Cone(u!, w±), L 2 = Cone(v 2 , ^2)1 we have the following equation of characteristic 
functions: 

Xk x + XKv = XLi + XL 2 i 

and hence 

M (ifi)+M(^ V ) = M(£i) + M(£2). 
Now the desired result follows from the fact that /i(L)(0) = 1/4 for orthogonal 2-dimensional cones L. To 
see this note that L and its image rot T / 2 (^) under rotation by it/2 form a half-space, and rotation is an 
isomorphism of V that preserves the inner product, hence fj,(L)(0) = /i(rot w / 2 (i))(0). 

Now suppose ^ is induced by a complete flag in V* . This flag is determined by a line U in V*, and ty* 
is induced by the dual flag, determined by the dual line U* in V. Now choose an identification i : V — > V* . 
Under this identification U pulls back to a line perpendicular to U* . It follows that 

^*(if v )(0) = A1 *(rot w/2 z- 1 (/f v ))(0). 

But rot 7r / 2 i _1 (_ftT v ) and K fit together to form a half-space. Hence 

M*(^)(0) + M *(rot 7r/2 r 1 (^ v ))(0) = i 

The result follows. □ 

Example 6.5. Consider the triangle P from Example 15 .61 with vertices v$ = (0, 0), V\ = (1, 0), w 2 = (0, 1). For 
^ corresponding to the standard inner product, recall that the constant terms of the y, of the vertex cones 
are 1/2, 3/8, and 3/8. According to the discussion above, we have that the constant terms of the is's are 
found by subtracting from 1/2: 

u(K )(0) = \ K#i)(0) = \ K# 2 )(0) = \. 

One can then check that the volume of this triangle is given by summing v(F)(0) values times the number 
lattice points in the relative interior of F. Indeed, since the only lattice points are the vertices, one has 

Vol(P) = - + - + - = -. 
V ; 4 8 8 2 
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For the triangle from Example l5.8[ with vertices = (0, 0), V\ — (2, 0), v 2 — (0, 1), one similarly computes 
the v(K) constant terms (again using the standard inner product): 

f(*o)(0) = \ K^i)(0) = ^ v{K 2 ){Q) = ~. 

This time there is one lattice point in the relative interior of an edge. The corresponding v constant term is 
1/2, according to Proposition 16.31 Thus, we check 
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Abstract. A local lattice point counting formula, and more generally a local Euler-Maclaurin formula 
follow by comparing two natural families of meromorphic functions on the dual of a rational vector space 
V, namely the family of exponential sums (S) and the family of exponential integrals (I) parametrized 
by the set of rational polytopes in V. The paper introduces the notion of an interpolator between these 
two families of meromorphic functions. We prove that every rigid complement map in V gives rise to an 
effectively computable Si-interpolator (and a local Euler-MacLaurin formula), an IS-interpolator (and a 
reverse local Euler-MacLaurin formula) and an IS°-interpolator (which interpolates between integrals and 
sums over interior lattice points.) Rigid complement maps can be constructed by choosing an inner product 
on V or by choosing a complete flag in V . The corresponding interpolators generalize and unify the work 
of Berline-Vergne, Pommersheim-Thomas, and Morelli. 
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1. Introduction 



1.1. What is a local lattice-point counting formula? The relationship between sums and integrals 
has been of interest to mathematicians since the ancient Greeks. The classical Euler-Maclaurin formula, 
discovered in the first half of the eighteenth century shortly after the development of modern calculus, may 
be viewed as a relationship between the sum of a function over the lattice points in a one-dimensional 
polytope with integer vertices and the integral of the function over the polytope. One naturally asks the 
same question in higher dimensions: Given a polytope P in an rt-dimensional space V equipped with an 
7i-dimensional lattice A and a function / on V, can one express the sum of / over the lattice points in P in 
terms of the integral / over P? In such a formula, one would expect a main term involving the integral of 
/ over P as well as correction terms involving the integrals of / over the proper faces of F C P. 

In the simplest case, let us suppose that / is a constant function. Then, the question becomes that of 
expressing the number of lattice points in P in terms of volume of P and the volumes of the faces F of P. 
If P is 2-dimensional, the celebrated Pick 's formula [Pi] 



expresses the number #(P) of lattice points in a convex lattice polygon in terms of its area A = vol(P), and 
the number of lattice points b — vol(9P) of its boundary. For example, we have: 



Here and throughout, we follow the usual convention that all volumes vol(P) of faces are normalized so that 
a lattice basis of AflL, where L is the linear space parallel to F, has volume 1. Unfortunately, Pick's formula 
is not local, due to the presence of the term 1. Here, locality means that for each face F of P, the formula 
contains a term that is the volume vol(P) multiplied by a coefficient that depends only on the supporting 
cone Supp(P, F) to P along F. The supporting cone is defined as the union of rays whose endpoint is in 
F and which remain in P for a positive distance. See Section [2~T1 below for a precise definition. McMullen 
[McMj proved the existence of local lattice point counting formulas. More precisely, he proved the existence 
of a function /i from rational, convex cones to rational numbers such that for any integral polytope P, the 
number $=(P) of lattice points in P is given by 



Here the sum is taken over all faces F of P, and vol(P) denotes the volume of the face F. 

Part of the difficulty constructing and computing a function fi that satisfies ([1]) is that /i is far from unique. 
The second author and Thomas gave an explicit construction of a rational valued function \x satisfying ((T|) , 
given a fixed complement map, a notion introduced in Thomas's thesis (cf. |Thj ) . A complement map 
is a systematic choice of complements of linear subspaces of a vector space; see Section 12.41 for a precise 
definition. Also note that all of the complement maps in this paper will be rigid. Two natural ways to get a 
rigid complement map are to choose: (1) an inner product, or (2) a complete flag. Given an inner product 
or a complete flag, |PTj construct effectively a map fi that satisfies Equation (fTJ). For the triangle depicted 
above, the Pommersheim- Thomas values of fi at the supporting cones to the vertices for an inner products 



#(P) = A+-b + l 




#(P) =4, A = 1,6 = 4. 



(1) 




F 
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or complete flags are as follows: 
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d\ -4^2+2^ 
6d 2 (d 2 -2d 1 ) 




1 
4 



9 dl+3d!d 2 +dl 



Udl-7d 1 d 2 +dl 
12di(2di-<J 2 ) 



20 12did 2 



The values on the left arise from the standard inner product on Z 2 , according to |PT( Cor.l], and may be 
computed by multiplying out the Todd polynomial in the ring presentation given in [PTl Prop.2]. Readers 
wishing more details on this should consult Section 15.21 The values in the figure on the right are those 
arising from the choice of the complete flag in V* = 1? whose i-dimensional subspace is spanned by the 
point (di,d,2). These values can be computed using the method outlined in |Mo[ p. 198]. They can also be 
computed by multiplying out the Todd polynomials using |PT1 Thru. 3]. One feature of this construction is 
that if one chooses a complement map arising from a complete flag, one recovers exactly the lattice point 
formulas of Morelli |Mo] , who gave a function /i that satisfies Equation (TT|) and takes values on the field of 
rational functions on a Grassmannian. For example, in the triangle on the right above, the values shown 
are rational functions on the Grassmannian of 1-dimensional subspaces of V*. The construction of |PT] is 
based on the theory of toric varieties, and gives an answer to a question of Danilov about the existence of a 
local expression for the Todd class of a toric variety. 

1.2. What is a local Euler-MacLaurin summation formula? Returning to the Euler-Maclaurin ques- 
tion, Berline and Vergne constructed in |BV] an explicit local Euler-MacLaurin formula for the sum of a 
polynomial function / over the lattice points P D A of an n-dimensional polytope P in a rational vector space 
V with lattice A. The Berline- Vergne formula has the form 



where the sum is over the set of faces F of P, and D(P, F) an infinite-order constant-coefficient differential 
operator D(P,F) that depends only on the supporting cone Supp(P, F). Equation ([2j is a generalization of 
(TT]). Indeed, if D(P,F) satisfies (|2|), and we define ^(Supp(P, F)) to be the constant term of D(P,F), then 
the local lattice point formula (JTJ holds. 

As in McMullen's case, the infinite order differential operators D(P,F) are not uniquely determined by 
([21). The construction of Berline- Vergne requires an inner product on the vector space V, and their results 
apply to rational polytopes, rather than just integral polytopes. Their construction associates to each cone 
K in V a meromorphic function (J>(K) on the dual space that is regular at the origin. For cones of dimension 
at most 2, the value of this function at (the constant term of the operator D(P,F)) coincides with the fi 
constructed in [PTj for the special case of complement maps arising from an inner product. In this way, for 
the example triangle above, the Berline- Vergne construction recovers the the numbers j, A arising out 
of the inner product case of the Pommersheim-Thomas construction. This fact that this coincidence holds 
is new to this paper; see Theorem [I] properties (5) and (6). 

1.3. An informal presentation of the results of this paper. The local Euler-MacLaurin formula @ is 
a consequence of a relationship between the integral and the sum of an exponential function over a polytope. 
In this paper, we introduce the concept of an interpolator between the families of exponential sums (S) and 
exponential integrals (J) over rational polytopes in a rational vector space V. We now discuss our results 
about interpolators informally, leaving the precise definitions and statements for Sections [2] and [3l 

If P is a rational polytope, or more generally a rational polyhedron (region, not necessarily compact, 
defined by linear inequalities with rational coefficients) in a rational vector space V, one can associate to 
P two important meromorphic functions, the exponential sum S(P) 6 Ai(V*) and the exponential integral 
I(P) £ Ai(V*) where A4(V*) is the algebra of meromorphic functions on the dual space V* <g>C. The values 



(2) 
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of these functions at a point £ e V* ® C are given by: 

(3) S(P)(£)= £ e«'*\ I{P){Q=fe^dx 

xePnA Jp 

provided \e^' x ^ | is summable (resp. integrable) over P. Here, the integral is taken with respect to the relative 
Lebesgue measure on aff(P), normalized so that a basis of A n W, where W is the linear subspace parallel 
to the affine span of P, has volume 1. The fact that Equations (|3J) define meromorphic functions, as well as 
the precise characterization and properties of the functions S and /, is essentially the content of Lawrence's 
theorem |Laj , reviewed in Section 12.21 below. 

The local Euler-Macluarin formula ([2]) follows in a straighforward manner from a formula of the following 
shape for a rational polyhedron P 

(4) S(P)=5>(Supp(P,F))J(F) 

F 

where [i is a function on the set of cones in V with values in A4(V*). We call such a function an SI- 
interpolator, or simply an interpolator. The main result of this paper (TheoremQJ states that a complement 
map on the vector space V gives rise in a natural way to an effectively computable Si-interpolator on V, 
and hence a local Euler-Maclaurin formula of the form [5] see Theorem [5J 

Thus, in particular, an inner product on V or a complete flag in V gives rise to a local Euler-Maclaurin 
formula. These interpolators have interesting connections with previous results. For complement maps 
arising from the choice of an inner product, one recovers the /i constructed in BV . In addition, we show 
that the values of these functions fi at are given by the functions \x in |PT] for cones of dimension at most 
2, and conjecture that these coincidences hold in all dimensions. In the case of complete flags, one obtains 
interpolators /i that are entirely new, as can be seen by calculations in dimension 2. In this case, the values 
of the constant term fj,(K)(0) coincide with values constructed by Morelli for cones of dimension at most 2, 
and conjecturally in all dimensions. One can also vary the chosen flag, thus associating to each cone K a 
meromorphic function fi(K) on F1(V*) x V* , where F1(V*) is the complete flag variety of V. This function 
fi(K) is naturally defined and effectively computable independent of any choices; see Theorem[3l In Morclli's 
work, he constructs what amounts to the constant term /i(if)(0) of this meromorphic function. Working 
dimension-by-dimension, Morelli chooses to view this term as a rational function on a Grassmannian, rather 
than on the entire flag variety. 

In Section 13. 2\ we show that a complement map on V also leads naturally to the construction of a local 
IS-interpolator, expressing the exponential integral I(P) in terms of the exponential sums S(F), for F a 
face of P; see Theorem 2J This construction, which works only in the case of integral polyhedra, allows one 
to obtain a reverse local Euler-MacLaurin formula expressing the integral of a polynomial function over a 
polytope in terms of the sums of the function over lattice points in the various faces of the polytope; see 
Theorem [6l 

Finally, we show how a complement map yields an IS°-interpolator, where we use S° to denote the sum 
over interior lattice points of a polytope. Part of the interest in IS°-interpolators lies in an observation of 
Morelli, who asked essentially (in the language of the present paper) if the same function can simultaneously 
serve as the constant term of an Si-interpolator and as the constant term of an IS°-interpolator on the dual 
space. He observed that the constant terms that he constructs do exactly this in dimensions at most 4. We 
give a proof of these coincidences in dimensions 1 and 2, and extend these results to the case of complement 
maps arising from inner products. Our arguments are essentially geometric, in contrast to Morelli's rather 
algebraic argument. (See Sections 13.31 and l6l for a complete discussion.) 

Section [5] contains explicit computations of the functions (J,(K) for some cones K of dimensions at most 
2. We compute these functions for complement maps arising both from inner products and from complete 
flags. In Section 15721 we exhibit Equation@]in the case of two triangles including the triangle depicted above, 
for both the inner product and complete flag cases, matching the constant terms with those constructed by 
Pommershcim- Thomas and Morelli shown in the above figure. 
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2. POLYTOPES, EXPONENTIAL SUMS / INTEGRALS AND INTERPOLATORS 

2.1. Polytopes. We will use standard terminology for rational polytopes, polyhedra, cones and faces, fol- 
lowing for example |Fu| . We will fix a finite dimensional lattice A, isomorphic to Z fe for some k £ N, and 
consider the rational vector space V = A ®z Q- All polyhedra will be rational and convex, and all cones will 
be rational, polyhedral and may be affine. A polytope is a compact polyhedron. If P is a polyhedron and F 
is a face of P, then the tangent cone Tan(P, F) is defined by picking a point x in the relative interior of F 
and looking at all the directions that one can go and stay in P: 

Tan(P, F) = {v G V\x + ev e P for small e > 0} 

This cone is independent of the choice of x and contains the origin. Let aff(P) denote the affine span of F, 
i.e., the smallest affine subspace of V that contains F. Let lin(P) denote the linear subspace of V parallel to 
aff(P). lin(P) is the maximal linear subspace of V contained in Tan(P, F). The supporting cone Supp(P, F) 
is the tangent cone translated back to its original position: 

Supp(P, F) = Tan(P, F) + x 

for x in the relative interior of F. 



2.2. Exponential sums and integrals. Let V* = Hoitiq(V,(Q)) denote the dual vector space to V, and 
(•, •) : V* x V — > Q denote the natural evalutation pairing on V. By a meromorphic function on V* we 
will mean a meromorphic function on the complexified dual space V£, where Vc = V <E)q C denotes the 
complexification of V. Let M.(V*) denote the algebra of meromorphic functions on V* and Ai r (V*) denote 
the subalgebra of meromorphic functions regular (i.e., analytic) at the origin. Note that if W is a rational 
subspace of V, then A' = A n W is a full rank lattice in W which naturally induces a relative Lebesgue 
measure on W, normalized so that a basis of A' has volume 1. Likewise, every rational affine subspace of V 
has a natural measure induced by the lattice. 

Given a polyhedron P in a rational vector space V with lattice A as above, there are two important 
meromorphic functions S(P) and I(P) with rational coefficients in V* , namely the generating function of 
exponential sums and exponential integrals. 

These functions are uniquely characterized by the following properties. 

(Al) If P contains a straight line, then S{P) = I(P) = 0. 

(A2) S (resp. I) is a valuation (resp. a simple (or solid) valuation). That is, if the characteristic functions 
of a family of polyhedra satisfy a relation J2i r iX(Pi) = 0> then the functions S(Pi) satisfy the 
relation r iS(Pi) = (resp. restrict the sum to those Pj that do not lie in a proper affine subspace 
ofV.) 

(A3) For every s € V and every £ £ V* we have 

(5) I(s + P)(0 = e^I(P)(0, 
and 

(6) S(s + P)(0=e^S(P)($, seA. 

(A4) If £ S V* is such that |e^ ,:c '| is integrable (resp. absolutely summable) over P, then 

(7) J(P)(0 = f e^dm P (x), S(P)(0 = £ e «'»> 

Jp xePnA 



Supp(P,F) 
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where dmp denotes the relative Lebesgue measure on aff(P). 

The existence of the functions / and S is by no means obvious, and is essentially the content of Lawrence's 
theorem; [La] . (Al)-(A4) determine the value of I and S on a simplicial cone as follows. 

(A5) If K = Cone(wi, . . . , v/.) is a simplicial cone generated by k independent vectors v%,...,Vk, (where 
k = dim(V) in the case of S) and □ («!, . . . , Ufc) denote the parallelepiped generated by vi,...,Vk, 
we have: 

is) w(o-(-i)' w ^|---.'* )) . w(o-f £ Hn^i^ 

lli=l\s» w »/ \xeD(Di,...,'K fc )nA / i=l 

Since any polyhedron can be subdivided (or virtually subdivided, using a linear combination of characteristic 
functions) into a union of simplicial cones, together with Equation (jSJ) above, one obtains an algorithm to 
compute I(P) and S(P). The complexity of such an algorithm is discussed in [Bvj . The next theorem which 
computes I(P) and S(P) of a polyhedron in terms of the tangent cones of its vertices was obtained by Brion 
(using toric varieties) and Lawrence (using combinatorics) independently. 

Proposition 2.1. [Brl ILa] Let P be a polyhedron in V. Then, 

(9) I(P)= Y, J(Supp(P,«)), S(P)= Y S(Snpp(P,v)) 

vev(P) veV(P) 

where the sums are over the set V(P) of vertices of P. 

2.3. Si-Interpolators. So far, our discussion of the meromorphic functions S(P) and I(P) was parallel but 
independent: our formulas did not mix 5* and /. It is natural to ask for an interpolation between S(P) and 
I(P). Indeed, if P is a polytope, then letting £ = reduces to a classic problem: expressing the number of 
lattice points in P in terms of the volumes of the faces of P. Let C(V) denote the set of cones of V. 

Definition 2.2. An Si-interpolator (or simply, interpolator) on V is map 

H-.C(V) -^M{V*) 

such that for any rational polyhedron P in V we have: 

(10) S(P)= Y i»(Supp(P,F))J(F) 

Fer(P) 

where the sum is over the set T(P) of all faces of P. An interpolator is regular if fj, takes values in M. r (V*). 

For the next definition, recall that if W is a quotient of V with projection tt : V — > W, then there is a 
natural map C(W) — > C(V) given by K M> 7r _1 (AT). Moreover, W* is naturally a subspace of V* , and there 
is a restriction map M(V*) — ► M(W*). 

Definition 2.3. An Si-interpolator /i on V is hereditary if for every rational quotient W of V, the composition 
ft, : C(W) — > M.(W*) given by the following diagram 

C(W) C(V) 



M(W*) M(V*) 

is an interpolator on W. 

The next lemma reduces the checking of the interpolator equation (|10[) to the case of cones. 
Lemma 2.4. /i is an interpolator if and only if it satisfies (|10|) for all cones P in V . 
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Proof. This follows from Proposition 12. II Indeed, we have: 

S(P) = ]T S(Supp(P,«)) 
vev(P) 

= E E M(Supp(Supp(P,«),F / ))/(F / ). 

veV(P) F'GJ r (Supp(P,i))) 

Since every face of Supp(Supp(P, v),F') is a face of Supp(P, P) for some face P of P, it follows that 

S(P)= £ M(Supp(P,P)) ^ J(Supp(F,t>))= ^ /i(Supp(P,P))7(P). 
FeJ r (p) «eV(F) fgf(p) 

The result follows. □ 

2.4. Rigid Complement maps. A rigid complement map gives us a systematic way to extend a function 
on a linear subspace of a vector space to the entire vector space. This notion was introduced in the thesis 
of Thomas |Thj and used in |PTj . where it is shown that a choice of complement map leads naturally to 
a local formula for the number of lattice points in a polytope and a local formula for the Todd class of a 
toric variety. In our paper, we will only use the notion of a rigid complement map. The curious reader may 
consult |PT[ ITh] for the definition of the general complement maps. 

Definition 2.5. A rigid complement map on V* is a map \& from a collection £* of linear subspaces in V* 
to the set of linear subspaces of V* satisfying two properties: 

(a) For every U G £*, *([/) is complementary to U; i.e., n U = {0} and *(P) + U = V* . 

(b) If E7i C J7 2 C V* with C/i, I7 a £ £*, then *({/ 2 ) c *(E7i). 

There are two easy ways to construct rigid complement maps: (a) an inner product Q on V* , and (b) a 
complete flag L on V*. 

Lemma 2.6. (a) An inner product Q on V* defines a complement map with domain £'3 the set of all 
subspaces of V* such that 

#(E7) = CP 1 

where [7 1 - denotes the subspace of perpendicular to U C V* under the inner product Q. 

(b) A complete flag L = (Lq, . . . , L n ) on V* satisfying Lq = {0} C L\ C • ■ • C L n — V* defines a complement 

map on V* by 

<3/ L (P) = i-„-dim([7) 

for any linear subspace U of V* which is complementary to £-n-dim((/)- This complement map is only defined 
generically, i.e., for subspaces that meet the flag generically. 

Our main theorem, Theorem [1] below, requires cones that are generic with respect to a rigid complement 
map $ . Let us define what generic means. Fix a rigid complement map ^ on V*. 

Definition 2.7. (a) A quotient W of V is ^-generic if W* G £*, where W* C V*. 

(b) A subspace U of V is ^-generic if V/U is 'J'-generic. 

(c) A cone K of V is \& -generic if lin(P) is Vt-generic for every face P of -ftT. 

(d) A polyhedron P of V is ^ -generic if for every face P of P, Supp(P, P) is ^-generic. 

The next lemma states that the notion of a complement map on V* is hereditary, i.e., a complement map 
on V* gives rise to a complement map on every W* , where W is a ^-generic rational quotient of V. 

Lemma 2.8. If is a rigid complement map on V* and W is a ^-generic rational quotient of V , then 
there is a natural complement map W on W* defined by 

#(Z7) = $((/)nr 

for all U C W* such that ^(U) (1 is complementary to U . 

The next lemma gives the promised extension of functions on a linear subspace W* to ones on the entire 
space V*. 
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Lemma 2.9. I/vp is a rigid complement map on V* and W is a -generic quotient ofV, then the decom- 
position 

V* = W* ®^{W*) 

defines a unique linear projection map 

7T : V* — > W* 

that annihilates ^(W*) and is the identity on W* . This allows any function on W* to be extended to V* . 

We now come to a main definition. Fix a rigid complement map $ on F. Let '(V) denote the set of 
^-generic cones of V. 

Definition 2.10. A ^-compatible interpolator fi on V is a function \i : C*(V) — > A4(V*) satisfying 

(a) Equation (fTU| for all VP-generic polytopes P, 

(b) \i is ^-hereditary, i.e., for every ^-generic quotient W of V, the composition p : C(W) — > M(W*) 
given by the following diagram 

C*(W) C*{V) 

(11) a m 

M(W*) - M(V*) 

satishes Equation (|10p for all ^-generic polytopes P in W . 

(c) Moreover, the following diagram commutes: 

C*(W) C*(V) 

(12) a » 

M{W*) —2— M(V*) 

where tt is given in Lemma 12.91 

3. Statement of the results 

3.1. SI- interpolators and local Euler-MacLaurin formula. Our main result is that a rigid complement map 
& on V* determines algorithmically a unique ^-compatible interpolator /i* on V. This generalizes results 
of [BV| IMo| IPTj . The proof uses the techniques of |BVj . which were a motivation and inspiration for us. 

Theorem 1. (a) If V is a rational vector space and "J is a rigid complement map on V* , there is a unique 
^/-compatible interpolator fi 9 , 

(b) In addition, fi satisfies the following properties on the collection of Vf-generic cones on V: 

(1) (Additivity) If the characteristic functions of a finite collection of cones Ki with vertex v 6 V satisfy 
the relation ^ rix(Ki) = 0, then the functions fj, (Ki) satisfy the relation J^. rj/i • (-Ki) = . 

(2) (Lattice Invariance) If v € A, then ^ (v + K) = ^ (K). 

(3) (Isometry Equivariance) If g is a lattice-preserving linear automorphism of V and g* is the inverse 
transpose, then (g(K))(g*($) = M*W(0- 

(4) (Regularity) ^* is regular at £ = 0. 

Moreover, /z* generalizes previous results of [Moj and [BVj : 

(5) (Constant Term) Let K be & top dimensional cone in V with vertex at zero, and suppose that V 
has dimension at most 2. Then for arising from an inner product or a complete flag, the constant 
term of ^ (K) agrees with the Todd class coefficient /i(AT v ) of the dual cone K v that depends on 
\& , constructed in [PTi Cor.l]. In particular, for VE' coming from a complete flag, the constant term 
of /z*(-K) gives Morelli's formula |Mo) . 
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(6) (Inner Product) For W coming from an inner product Q* on V* , the function /i* ' (K) agrees with 
[j.(K) constructed in |BVj . where K is the image of K in V/V(K), where V(K) denotes the linear 
subspace parallel to the largest affine subspace contained in K. 

There are two corollaries of Theorem [TJ a local version of the Euler-MacLaurin formula for polytopes 
(promised in the introduction), and a meromorphic function on F1(T^*) x V* associated to a pointed cone K 
in V, where F1(V*) is the variety of complete flags of V. 

To formulate these results, fix a rigid complement map ^ on V , a ^-generic rational convex polyhedron 
P and a face F of P. Consider the Taylor series expansion of /i (Supp(P, P))(£) and the corresponding 
differential operator (of infinite order, with constant coefficients): 

(13) D*(P,P) = /i*(Supp(P,P))(9 x ). 

Theorem 2. Let ^ be a complement map on V* , let P be a ^ -generic rational polytope in V , and let h(x) 
be a polynomial function on V . With the above notation, we have 

(14) h ^= E / D*(P,F)-h. 
In particular, if h(x) — 1 then 

(15) #(P)= M*(Supp(P,P))(0)vol(F) 

FEF(P) 

where #(P) (resp. vo\(P)) denotes the number of lattice points (resp. the volume) of P. 

For concrete examples, see Section [5] Equation (|15j) computes the number of lattice points of a polytope 
as a weighted sum of the volume of its faces. 

Suppose now that ^ comes from a flag L. Letting L vary, we obtain the following. 

Theorem 3. Given a pointed cone K in V there exists an effectively computable meromorphic function 
fJ,(K) in F1(V*) x V* which is regular in a Zariski open subset o/Fl(V*) X {0}. 

For an algorithmic computation of n(K) for cones K of dimension at most 2, see Section [5] 
The next conjecture identifies the constant term of a ^/-compatible interpolator with the one of |PT| . We 
post it as a conjecture for now, and hope to discuss it in a future publication. 

Conjecture 1. Part (5) of Theorem{l\ holds in all dimensions. 

Question 1. Is there any relation between the rational functions fi(K) in F1(V*) X V* and the quantum 
cohomology ring on Fl(l^*)? 

3.2. IS-interpolators and reverse Euler-MacLaurin formula. The definition of the Si-interpolator leads in a 
natural way to the notion of an IS-interpolator, which may be a useful notion to numerical approximations 
of integrals by sums. The next notion requires us to restrict attention to lattice polyhedra, i.e., polyhedra 
whose vertices are points of the lattice A. Below, a lattice cone K in V is a rational cone such that UDA ^ 0, 
where U is the largest affine subspace contained in K . 

Definition 3.1. An IS-interpolator is map A, from the set of all lattice cones in V to Ai(V*) such that for 
any lattice polyhedron P in V we have: 

(16) J(P)= A(Su PP (P,P))S(F) 

where the sum is over all faces F of P. An IS-interpolator is regular if A takes values in Ai r (V*). 

Hereditary and ^-compatible IS-interpolators are defined in an analogous way to Definition 12.31 and 12. 101 
The statement and proof of Theorem [1] holds with minor modification for IS-interpolators. 
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Theorem 4. (a) If V is a rational vector space and \1/ is a rigid complement map on V* , there is a unique 
^-compatible IS-interpolator A*. 

(b) In addition, satisfies properties (2)-(5) of Theorem [T] and the following version of additivity: 

(1) (Additivity) If if is a cone that does not contain a linear subspace and is subdivided into a finite 
union of cones if*, then A* (if) = Ei:dim(if»)=dim(K) >?{ K i)- 

One can obtain directly a formula for A* in terms of the interpolator /x* of Theorem [1] as follows, by 
observing that /i*(Supp(if, if )) = 1. 

Theorem 5. //if is a positive (resp. zero) dimensional lattice cone then 

(17) Yl A*(iV*(Supp(K,F)) = 0, f^.J A*(if) = 1. 

The interpolator A* gives a reverse Euler-MacLaurin summation formula. Fix a rigid complement map 
^ on V, a ^-generic lattice polyhedron P and a face F of P. Consider the Taylor series expansion of 
A*(Supp(P, P))(£) and the corresponding differential operator (of infinite order, with constant coefficients): 



(18) A*(P,P) = A*(Supp(P,P))(cy 

Theorem 6. With the above assumptions, for every lattice polytope P and every polynomial function h(x) 
on V we have 

(19) / h(x)= ]T E {^(P,F)-h)(x). 

p FeJ r (p) i£fnA 

In particular, if h{x) = 1 then 

(20) vol(P)= J2 A*(Sup P (P,P))(0)#(P) 

where #(P) denotes the number of lattice points of F. 

Equation (j2"0|) computes the volume of a lattice polytope as a weighted sum of the number of lattice points 
of its faces. A formula of this type was first written down by Morelli, |Mo[ Eqn.(5)]. 

3.3. IS°-interpolators and Morelli's work. When ^ comes from a complete flag, Equations (fT5|) and 
(|20|) arc similar with some results |Mo| Eqn.(5),Eqn.(6)] of Morelli. To explain this, let us introduce the 
variant S° of the exponential sum function defined by: 



(21) S°(P) = S{P°) 

where P° denotes the (relative) interior of a polyhedron P. The set P° is not a polyhedron itself, however 
it is a virtual sum of polyhedra. Then, we can talk about IS°-interpolators v on lattice cones. Theorems 2] 
and [5] have the following analogue. 

Theorem 7. (a) If V is a rational vector space and \& is a rigid complement map on V* , there is a unique 
^-compatible IS -interpolator V s . In addition, 

(b) (Additivity) If if is a cone that does not contain a linear subspace and is subdivided into a finite union 
of cones K l7 then v*(K) = £ rdim(/f! )=dim(if) v 9 (Ki). 

(c) v satisfies properties (2)-(5) of Theorem[T] 

(d) If if is a lattice cone then 

(22) J2 ^(^)M*(Supp(if,P)) = l. 

FeF(K) 
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Fix a lattice polytope P and a face F of P. If \& comes from a complete flag L in V* as in Lemma \TM and 
we vary the flag, we obtain two rational functions i/(Supp(P, F)) and /i(Supp(P, F)) on F1(V*) x V* which are 
regular at a Zariski open subset of F1(V*) x 0. Thus, we can consider the constant terms i/*(Supp(P, P))(0) 
and /x*(Supp(P, P))(0) which are rational functions on Fl(V*). 

On the other hand, Morelli constructs a pair of rational functions ^(tipP) and /j,k{$FP) (where k is 
the dimension of lin(Supp(P, F))) that appear respectively in Equations (5) and (6) of |Moj (unfortunately, 
Morelli denotes the two functions with the same notation, although he clearly recognizes that they are 
distinct functions). Morelli's construction uses K-theory, localization and the Bott residue theorem. 

The following proposition is a re-expression of the flag case of Theorem [1] part (5). 

Proposition 3.2. With the above notations, for polytopes P of dimension at most 2, we have: 
(23) KSu PP (P,P))(0) = M$fP), /i(Su PP (P,P))(0) = nuippP). 

In [Mol p. 191] and also in [Mol Thm.5], Morelli observes that 



(24) »(K)(0) = K^ V )(0) 

for lattice cones K of dimension 4 or less, and that Equation ([24| fails in dimensions more than 4. This 
mysterious and seemingly deep coincidence is easy to explain from our point of view, for cones of dimension 
at most 2. The rational functions v(K ) and n(K v ) are distinct, in fact they are functions on distinct spaces 
(Fl(y*) x V* and Fl(y) x V, respectively). Their constant terms coincidentally agree in small dimensions. 
For a discussion, see Section 16.21 

4. Proofs 

4.1. Proof of Theorem [TJ Fix a rigid complement map $ on a rational vector space V with lattice A, and 
consider a ^-compatible interpolator /i* and a Vf-generic rational polytope P. We will show by induction 
on the dimension of P the existence and uniqueness of /i*(P). All cones in this section will be assumed 
to be ^-generic. By Lemma 12.41 it suffices to assume that P = K is a cone in V. To avoid confusion, we 
denote the map /1 : (V) — > A4(V*) by [iy and for every 'I'-generic quotient W of V, we denote the map 
fd : C*(1F) — ► M{W*) by ft w . 

If K is pointed with vertex v, then we can single out the contribution from the 0-dimensional face of K 
in Equation (|10|) , and together with hereditary property of Equation ([TT]) , it follows that 

(25) $(K)(t) = e-^l ]T e«'*>- £ 4 /lhl(F) (Sn PP (K,F))(n(0) f e^dm F (x) 

\x£KnA _F,dim(_F)>0 F 



where Supp(i^, F) denotes the image of Supp(P, F) in the quotient space V/lm(F), and ir : V* — > 
(Vyiin(P))* is the projection using the decomposition V* = (F/lin(P))* © *((V/lin(P))*) given by the 
rigid complement map ^/ 

On the other hand, if if is a non-pointed cone in V, let V(K) denote the linear subspace of V parallel to 
the largest afhne subspace contained in K, and consider the image K of K in V/V(K), which is a pointed 
cone of dimension strictly less than the dimension of K . Equations (111|) and (|12p imply that 

(26) n v (K)(Q=4 /v{K) (K)(0 

where £ € V* , tt : V* — > (V/V(K))* is the projection using the decomposition V* = (V/V(K))* © 
^((V/V(K))*) given by the rigid complement map \& and £ £ (V/V(K))* is the image of £ under the above 
projection. 

Equations (f2l)|) and (f2l)]) uniquely define fj, v from for dim(Vy) < dim(V r ). On the other hand, when 
V is a 0-dimcnsional space, we have My({0}) = 1- This uniquely determines n v . An explicit computation 
of ^(K) for cones K of dimension 0, 1 and 2 is given in Section [5j 
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Additivity of /i* follows from the above inductive definition of /x* and the additivity (property (A2)) of 
the exponential sum and exponential integrals. This is presented in detail in the proof of BV, Prop. 15]. 

Remark 4.1. In fact, Berline-Vergne in |BV] work with the transverse cone T(K, F) of a face F of K defined 
to be the image of Supp(2f, F) under the projection map V — > V/\in(F). The transverse cone T(K, F) is 
always a pointed cone. In fact this gives a 1-1 correspondence 



(27) {cones in V} O {pointed cones in quotients of V} 

Using this correspondence, and an inner product Q on V, Berline-Vergne inductively construct /uP defined 
on the set of pointed cones on the quotients of V. For visual reasons, we prefer to work with the supporting 
cones Supp(if, F) rather than the transverse cones T(K, F). 

Going back to the proof of Theorem [TJ the lattice invariance and the isometry equivariance of //* follows 
easily by induction and the invariance of the I and S functions. 

The regularity of the meromorphic function /j, ' (K) at zero follows by induction and a residue calculation, 
discussed in detail in }BV1 Prop. 18], using the following fact about the exponential sum and integral functions. 

(A6) If if is a cone in V with primitive integral generators v±, . . . , Vk, then 

fe k 

(28) n&^WKO. l[{Z,Vi)I(K)(t) 

i=l i=l 

are regular functions and the residues of the meromorphic functions S{K) and I{K) along the 
hyperplane V\ = satisfy the equations 

(29) Res^ (S(K)) = -S(K), Res Vl {I{K)) = -I(K). 

(A6) follows in turn by additivity and Equation ((SJ which evaluates the I and S functions on simplicial 
cones. 

We defer the proof of Property (5) of Theorem [1] until the example section (Section [5]). See the proofs 
immediately following Examples 15.61 and 15.71 

Finally property (6) of Theorem Q] follows directly from |BV[ En.4] . 

4.2. Proof of Theorem [2j Fix a polynomial function h(x) on V and a polytope P on V. Consider the 
constant coefficient differential operator h(d^) which acts on the function £ i-^ e^' x ) by: 

(30) h(d 6 )e^ =h(x)e^' x l 
The definition ([7]) of S(P) and the above gives that 



.... ,^ ] ' x) 

xePnA 

Let ev denote the evaluation at £ = 0. It follows that 



(31) (evoh(d ( ))S(P) 
On the other hand, Theorem [1] gives that: 



xg-PnA 
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S(P)(0 = E /i*(T(P,F))(OJ(F)(0 
FeJ r (P) F 

= E / M*(T(P,F))(^)e«^^mp( a; ) 



Fgf(P) 

where the last equality follows from the definition (fT3|) of the differential operator D 9 (P, F). Applying the 
differential operator /i(c^) to both sides it follows that 

h(d £ )I(P)(Q = E / D*(P,F)h(x)e {s < x) dmp(x). 
FeT(P) Jf 

Evaluating at £ = 0, it follows that 

(32) (evoh(dt))S(P)= E / D*(P,F)-h. 

FeT(P) F 

Equations (|3 1 1) and (|32|) complete the proof of Theorem [2] 

4.3. Proof of Theorem|4l The proof of Theorem [T] applies verbatim, with the following observation, which 
explains the need for lattice polyhedra. 

If K is a rational pointed cone with vertex v for which (|16l) applies, then one of the terms in (fTo]) is 
F = {v}. In that case, T(K, {v}) — K and the corresponding term in (fTo) is ^(if )5({i;}). Now, (A3) and 
(A5) imply that 



e^' v ~> if v £ A 
otherwise. 



Thus, when K is not a lattice cone, then we cannot solve for v(K). In the case of SI- interpolators, as 
discussed in Section [2~4l the corresponding term of (fTOj) was fi(K)I({v}) and (A3) and (A5) imply that 

J(W)(0 = e«'»>. 

4.4. Proof of Theorems [5] and O Theorem [5] follows from the fact that the composition of an SI- 
intcrpolator with an IS-interpolator is an SS-interpolator, and the fact that an SS-intcrpolator is unique. 
More precisely, fix an Si-interpolator /i and an IS-interpolator A, and consider a lattice cone K. Then, on 
the one hand we have 

S(K)= E ^npp(K,F'))I(F'). 

On the other hand for every face F' of K, we have 

I(F')= E X(Supp(F',F))S(F) 

FeF{F') 

Substituting the second equation into the first, it follows that 

S(K)= E ^(Sup P (A',F'))A(Supp(F',^))5(F). 

F,F': FCF' 
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This motivates the following definition. Consider the function 

Supp(P, F) h+ (A o M )(Supp(P, F)) := ]T A*(Supp(if, F'))A(Supp(F', F)). 

F': FcF' 

It is easy to see that Xo/j, is an S'S'-interpolator. Moreover, if \x and A are ^/-compatible, so is their composition. 
On the other hand, there is a unique ^-compatible ^-interpolator. Thus (A o //)(Supp(P, F)) = Oresp. 1 
for F ^ P (resp. F = P). This concludes the proof of Theorem [5] 

Theorem [7] follows by an analogous computation, using the fact that the map (y o fi) is a ^-compatible 
SS°-intepolator, the uniqueness of such interpolators, and the fact that the map that sends every cone to 1 
is a ^/-compatible S , S' -intcpolator. 

5. Computations and Examples 

5.1. Computation of (jr for low dimensional cones. As was mentioned earlier, the meromorphic func- 
tions jit*(if) are effectively computable given a complement map \&. In this section we illustrate this, by 
explicitly computing for cones K of dimension at most 2. 

Proposition 5.1. Suppose that K = {0} in the vector space V = {0}. Then fi(K) is the constant function 
1, independent o/W. 

Proof. S(K) = I(K) = 1 forces fj,(K) = 1. □ 

If {vi, . . . , Vk} is a collection of vectors in V, let (vi, . . . ,Vk) = Q + «i + ■ ■ ■ + Q + Vk denote the cone spanned 
by those vectors, where Q + is the set of non-negative rational numbers. 

Proposition 5.2. Now suppose that V is one- dimensional and K — Cone(u) is a ray in V generated by a 
primitive vector u £ A. Then, independent of ^ , we have 

where 

1 1 1 z z 3 z 5 



^ B ^ Z ' 1 - e z ' z 2 12 ' 720 30240 

is the generating series of the Bernoulli numbers. 

Proof. Use the property of interpolators 

S{K) = n(Supp(K, K))I{K) + fi(Supp(K, 0))7(0), 

/i(Supp(i\T, K)) = 1 by the previous proposition and compatibility under quotients. In addition, we have 
S(K)(€) = 1/(1 - e«) and = The result follows. □ 

For the following proposition, recall that if K is a subset of V, then we define the dual K y by 



(34) K v = {w G V* | (w, v) > 0, for all v G K}. 

Also, if A is a lattice in V, then A* = {w G U* u) G Z, for all v G A} is the dual lattice in V* . Below, 
we will denote by /i L (resp. /jP) the ^/-compatible Si-interpolator where ^ comes from a complete flag L on 
V (resp. an inner product Q on V), as in Lemma 12.61 

Proposition 5.3. Suppose V is two-dimensional, and K C V is a half-plane with boundary a line U through 
the origin. 
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Let p G A* be the primitive generator of the ray K v . Suppose the 'J is defined by a complete flag L in V* , 
and let c be a generator of L*, where L\ is the one- dimensional subspace of the flag L (the only nontrivial 
subspace in the flag.) Then 

' (Ley 

Now suppose ^ is defined by an inner product Q on V . Then 



p Q (i<m = b 



Q(Z,p) 
Q(p, p) 



Proof. We use the ^-compatibility together with our one-dimensional formula above. To do so, we must 
compute the projection from tt : V* ^ W* where W = V/U. This projection is computed using the 
decomposition V* — W* "i>*(W*) given by the complement map Suppose the 1-dimensional subspace 
^(W*) C V* is generated by d. Then for £ E V* , find scalars uji,u> 2 such that 

(35) £ = ujip + cj 2 d, 

Then from ^-compatibility and the one-dimensional formula, we obtain 

/i(A-)(0 = S(wi). 

For coming from a complete flag, we pair both sides of Equation (|3"5"|) with the generator c. Since (d, c) = 0, 
we find 

u ( ^' C> 

Wl - -. r, 

and the formula for p L (K) follows. For \1/ coming from the inner product Q, we note that Q(d, p) — 0. Thus 
Equation ([33|) yields 

Q(p,p)' 

which completes the proof. □ 

Using the defining property of interpolators, we immediately get the following proposition, an explicit 
expression for p of a nonsingular 2-dimensional cone. 

Proposition 5.4. Suppose V is two-dimensional and K = Cone(ui, v 2 ) where v\, v 2 form a basis of A C V. 
Let Fi = Cone(ui) and F 2 — Cone(t>2), and let p\, p 2 G A* denote the primitive normals to F± and F 2 . Then 
for ^ coming from a complete flag L in V* , we have 

r ,<n „l (k) = i i__ | i R (&±\ | 1 

( ' Ml ' (l-e«.»0)(l_ e «^>) (t,Vl){Z,V2) + (£,Vl) \(P1,C)) + \( P2 ,c) 

For "J coming from an inner product Q on V* , we have 

{ ) Ml ' (l-e«.«i>)(l-e«^>) <e,«i><e,«a> \Q(pi,Pi)J (t,t») \Q( P 2,P2) 

Remark 5.5. Replacing 1/(1 — e x ) with B(x) — 1/x, it follows that Equations ([55)1 and (|3"T|) can be written 
in the form: 

T, s ,s >s 1 ( ( (f,c) \ ls \ 1 / / (£,C) 



* ( *> " + - B <«-»J + (n<^> 1 - B( «'"' ,) 

. b({{ ,„ i)) b((^» + ^(b(^1)- B (( ? ,„ 2 ),) 
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5.2. Examples. In this section we explicitly compute some low dimensional examples and match them with 
the ones given by Morelli, Pommersheim-Thomas and Berline-Vergne. 

Example 5.6. Consider the triangle P with vertices vq — (0,0), V\ — (1,0), v 2 = (0,1) in the lattice A = 
1? C V = Q 2 , where V* is also identified with Z 2 . 




Let ^ be the complement map coming from the complete flag L = (Lq, Li,L 2 ) in V* such that L\ = 
Span{((ii, da)}- To compute fi L of the supporting cones, we use lattice invariance to translate them to the 
origin. In other words, consider the tangent cones K{ = Tan(P, Vi) for i = 0, 1, 2, to P at V{. Specifically we 
have 

K Q = Cone((l, 0), (0, 1)), K x = Cone((-l, 1), (-1, 0)), K 2 = Conc((0, -1), (1, -1)). 
Then for G V*, we have 

r, , 1 1 1 ( dx \ 1 / d 2 

m = (1 _^. 1 )(1 _ e - e . ) + ^ - - ^) + j^-A^r 



^{K 2 ){£) =————-— + — — —B[ £i — — -£ 2 + ^ 



(l-eCi-fo)(l-e-&) 6(6-6) 6 V d 2 J 6-6 V di-da 
These functions are analytic at £ = (0,0). For example, ii l (Kq)(£) has the Taylor expansion 
L / df+j|+3^2 \ /-IV /-IV / Bjf+j^ + Srffjf V 

720d?d 2 J^ 1 + V 720d|di J^ 2 + " ' 
Denoting the constant term of each fj, L (Ki) by /Uq (-Ki), we find 

L d 2 + d 2 +3did 2 

Mo (^o) = - 



Mo(^i) 

^0 



12gW 2 
5d 2 - 5did 2 + d| 

12di(di -efa) 
5dl - 5oW 2 + d 2 



12d 2 (d 2 -d 1 ) 

in exact agreement with Morelli's example |Mo[ p. 198-199]. These rational functions can also be computed 
by multiplying the Todd polynomials using (PT1 Theorem 3] In particular, these three rational functions 
sum to 1 in agreement with Equation (1151) . 

In addition, one checks by direct computation that 

J2^ L {^pp(P,F))I(F) = S(P) = 1 + e? 1 +e« 2 . 

F 

in agreement with the definition of an interpolator. 

Indeed, this example provides a general proof of the coincidence of our construction and the flag case of 
|PT) for cones of dimension at most 2, and hence the coincidence with Morelli's formulas. 

Proof of QJ part (5) for Bags: In dimension 1, the |PT] values are always equal to 4, independent of 
complement map. Using Proposition 15.21 we see that this coincides with the constant term of /x. For the 
two-dimensional case, we first use the additivity of |PTj and of our construction (Theorem 1, Property (1)) 
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to reduce to nonsingular 2-dimensional cones. By functoriality, we may assume that K is the cone K from 
example 15. 61 generated by the standard basis of Z 2 . In this case, we note that the values of p,q(Kq) given in 
Example 15.61 match the values in |Mo[ p. 198-199]. By |PT[ Corollary 2], these values also match the values 
in |PT) . This completes the proof. 



Example 5.7. We consider the same triangle P from Example l5.61 but now with a complement map ^ coming 
from an inner product Q on V*, given by the matrix 

b ci 



We can use Proposition 15.41 to compute ^(Ki) for the supporting cones Ki of P. We find 

^(*o)(0 = T , KTfi 77T -TT + T B (^ + + T B U 2 + 

(1 - e>i)(l - e* 2 ) Ci6 6 V a J ?i V c 

" ( * l)(0 = (i-e^-^)d-^O + Mi^a) " 6 fl l 6 + c^J + 6^iA 

Qfrs \m 1 1 1 p/^ ,b\ 1 p ^<l-^l-^2-^2 

This time, we find the constant terms are given by 

Sac — ab — be 



\2ac 

ab + 4ac + 106c + 2b 2 + 5c 2 

12(ac + 26c + c 2 ) 
5a 2 + 26 2 + lOafe + 4ac + be 



12(a 2 + 2a6 + ac) 

These agree with the values of p, constructed in |PT1 Cor.l] for coming from an inner product. In 
particular, the three rational functions above sum to 1. For the standard inner product on V* = Z 2 , 
(corresponding to a = c = 1, b = 0), we find 

For completeness and comparison, we detail how these numbers arise, in a seemingly very different way, 
out of the construction of Pommersheim-Thomas. We begin, according to |PT[ Cor.l(iv)] by considering the 
outer normal fan of P. The rays of this fan are pa = (1, 1), pi = ( — 1, 0), and p 2 — (0, —1). We must then 
multiply the second-degree Todd polynomial T = J2i<j DiDj + i J2i Df in the ring given in |PT1 Prop. 2]. 
In this ring, one finds D 2 = + \D D 2 , D\ = D a Di, and D 2 — D D 2 . Hence 

T = -D X D 2 + \d q D 2 + Id D u 
4 o o 

from which one reads off the values of po(Ko) = j, po(Ki) = |, and po(K 2 ) = |. In fact, more generally, 
let K = cone(t>i, v 2 ) be any nonsingular cone in V with vertex at 0, and let Q be an inner product on V* . 
Then an easy computation similar to the above can be used to compute the value of po of [PT| for this cone 
K. Indeed, if K* — cone(pi, p 2 ) is the dual cone, then using presentation of |PT| Proposition 2], one finds 
that the coefficient of D\D 2 in T is given by 

(38) 1 Q(fii,Pn) (I , 1 \ _ 1 , Q*(«i,«a) 



4 12 \Q(j>i,pi) Q(P2,P2)J 4 12 V<2*(«i>«i) Q*{v2,v 2 ), 

where Q* denotes the dual inner product on V . 

Indeed, this calculation proves the inner product case of Theorem [TJ part (5), as we now show. 
Proof of [TJ part (5) for inner products: As in the flag case, the |PT] values are always equal to -|, 
independent of complement map, matching the constant term of the functions in Proposition 15.21 For the 
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two-dimensional case, we again use additivity and functoriality to reduce to the case in which K is the cone 
Kq from example 15.61 We then note that the formula of Equation [35] for the |PT| construction matches the 
values of /Iq(Ko) given in Example 15.71 This completes the proof. 

A comparison of the functions n®{Ki) of Example 15.71 with the functions ^{Ki) of Example 15.61 reveals 
that they are indeed different, confirming the fact that interpolators that come from flags are different from 
those that come from an inner product. 

Let us repeat the previous example with a different triangle. 

Example 5.8. Consider the triangle P with vertices vq = (0, 0), V\ = (2, 0), v 2 = (0, 1). 




As above, let VP be the complement map coming from the complete flag L = (Lq, Li, L 2 ) in V* such 
that LI = Span{(di, d 2 )}. Denoting, as before, the tangent cones at the vertex Uj by K i7 we see that 
Kq = Cone((l, 0), (0, 1)), and K\ = Cone((— 1, 0), (— 2, 1)) are nonsingular cones, whose \i values may be 
computed from Proposition 15.41 On the other hand, the cone K 2 — Cone((0, —1), (2, —1)) is singular. We 
can compute n L (K 2 ) by subdividing K 2 into two cones 

K 2 = d U C 2) d = Cone((0,-l),(l,-l)), C 2 = Cone((l, -1), (1, -2)) 

using the ray p = Cone((l, —1)). The additivity of /i implies that 

p L {K 2 )=p L (C 1 )+p L (C 2 )-p L (p). 

The terms on the right hand side can be computed using Propositions 15.41 and 15.21 Alternatively, one 
can compute /j, l (K 2 ) directly from the defining formula for interpolators. Letting F\ — Cone((0, —1)) and 
F 2 = Cone ((2, —1)) denote the one dimensional faces of K 2 , we get 

p L (K 2 ) = S{K 2 ) - I{K 2 ) - fi L (Snpp(K 2 ,F 1 ))I(F 1 ) - M L (Su PP (^ 2 , F 2 ))I(F 2 ) 

1 + *-* , 2 , 1 B (^_d^\ + 1 



(l- e -6)(i_ e 2 ei -6) 6(26-6) -6 V d 2 sz J 26-6 V d 2 -2d, 
The constant terms of \x L (Ki) are Morelli's rational functions 

L d\ +Mid 2 +d\ 



12did 2 



L lldj - ld x d 2 + d\ 

^ [Kl >- Ud 1 (2d 1 -d 2 ) 

L dl-M 1 d 2 + 24 

= _ 6d2(2dl _ d2) 

agreeing with the numbers shown in figure in the introduction. 

Again, one can verify by direct computation that 2<=n MoC-^i) = 1 an d 

^/x i (Supp(P,F))J(F) = S(P) = l + e 6 + e 2Sl + e« 2 . 

F 

Example 5.9. We consider the same triangle P from Example l5.81 but now with a complement map ^ coming 
from an inner product Q on V* , given by the matrix 

As in Example 15 .81 the function p® of the singular cone K 2 may be computed either by subdividing, or by 
direct use of the defining property of interpolators. We find that /jP(K 2 ) is given by 

1 + , 2 , J_ B / _ b \ + 1 B ( -ofc - 266 - &6 - 2c6 N 



(l-e-«=)(l- e 26-e 3 ) 6(26-6) -6 V « / 26-6 V a + 46 + 4c 



SUM-INTEGRAL INTERPOLATORS AND THE EULER-MACLAURIN FORMULA FOR POLYTOPES 



19 



The constant terms, in agreement with the construction of [PTj . are 

o , T , s 3ac — ab — be 



\2ac 



Q ab + 5ac + 46 2 + 256c + 22c 2 



12(ac + 46c + 4c 2 ) 
Q _ 2a 2 + 8ab + lac + 2b 2 + 2bc 
Mo ( 2) " 6( a 2 + 4a6 + 46c) 
which add up to 1. For the standard inner product, these coefficients are 

Again for comparison, let us detail how these numbers arise out of the construction of Pommersheim- 
Thomas. To apply |PT[ Cor.l(iv)], we first consider the outer normal fan of P. The rays of this fan are 
Po = (1,2), pi = (—1,0), and P2 = (0,-1). A nonsingular subdivision is achieved by adding the ray 
P3 = (0, 1). As before, we multiply the Todd polynomial T — J2i<j DiDj + i Yli &l m the ring given in 
[PTl Proposition 2]. One calculates D 2 = \{2D Q D 2 - 2D D 3 ), D 2 = 0, B\ = 2D D 2 , and D 2 = -2D D 3 . 
Hence 

T = \d x D 2 + ^D Q D 2 + ^D a D 3 + jDiAs, 

Using additivity under subdivision |PT1 Cor.l(iii)], one finds pq(Kq) = j, p (Ki) — and po(K 2 ) = 
j_ , l j_ 

20 ' 4 10' 

6. Comparison with Morelli's work 

In this section we compute values of v and discuss the relations Morelli noticed between the constant 
terms of p and v. 

6.1. Computation of v for low dimensional cones. For cones of dimension 0, 1 and 2 it is easy to 
compute v (K) explicitly, just as we did for p, (K) in Section [5j The proofs in this case are similar. 

Proposition 6.1. Suppose that K — {0} in the vector space V = {0}. Then v{K) is the constant function 
1, independent of ^ . 

Proposition 6.2. Now suppose that V is one- dimensional and K — Cone(u) is a ray in V generated by a 
primitive vector v G A. Then, independent of ^ , we have 

Proposition 6.3. With the notation of Provosition [573[ we have 
and 

Proposition 6.4. With the notation of Proposition \5~^j\ we have 



<^«i)«,«a> 
and 



«.«0 V (Pi,c)J l-e«.«») V (P2,c)J (i_ e «.*i>)(l-e«>«»)) 



^ V Q(Pi,Pi)J l-e«.«»> V Q(P2,P2)J ' (l-e«-"i>)(l-e«^>) 
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6.2. Morelli's coincidence in dimensions 1 and 2. In this section we will prove Morelli's observation 
for lattice cones of dimension at most 2, using the computations of v^(K) for cones K of dimension 
at most 2 from Section 16.11 Using these propositions, we can prove the following relation between fi and v 
for cones of dimension at most 2. The following theorem, proved by Morelli in the flag case, shows that the 
coincidence of constant terms holds also in the case of inner products. 

Theorem 8. Let K be a cone with vertex at in V , and assume that the dimension of V is at most 2. Let 
^ be a complement map on V induced by an inner product on V* or a complete flag in V* and let be the 
corresponding complement map on V . Then we have 

Proof. For cones of dimension 0, the constant term of both p, and v is always equal to 1. For cones of 
dimension 1, the constant term of both /i and v is always 1/2. Now let K have dimension 2. From the 
inversion formula of Theorem it follows that 

Hence, the desired conclusion is equivalent to 

M^)(o) + M^ v )(o) = i 

Suppose that fi is induced by an inner product. The inner product gives us an identification of V with 
V* , allowing us to consider K and K v as both living in the same space V. If K = Cone(ui, V2), then 
K y = Cone(wi, w%) with m orthogonal to Wi, i = 1,2, under the inner product. Furthermore, considering 
the orthogonal cones L\ = Cone(u!, w±), L 2 = Cone(v 2 , ^2)1 we have the following equation of characteristic 
functions: 

Xk x + XKv = XLi + XL 2 i 

and hence 

M (ifi)+M(^ V ) = M(£i) + M(£2). 
Now the desired result follows from the fact that /i(L)(0) = 1/4 for orthogonal 2-dimensional cones L. To 
see this note that L and its image rot T / 2 (^) under rotation by it/2 form a half-space, and rotation is an 
isomorphism of V that preserves the inner product, hence fj,(L)(0) = /i(rot w / 2 (i))(0). 

Now suppose ^ is induced by a complete flag in V* . This flag is determined by a line U in V*, and ty* 
is induced by the dual flag, determined by the dual line U* in V. Now choose an identification i : V — > V* . 
Under this identification U pulls back to a line perpendicular to U* . It follows that 

^*(if v )(0) = A1 *(rot w/2 z- 1 (/f v ))(0). 

But rot 7r / 2 i _1 (_ftT v ) and K fit together to form a half-space. Hence 

M*(^)(0) + M *(rot 7r/2 r 1 (^ v ))(0) = i 

The result follows. □ 

Example 6.5. Consider the triangle P from Example 15 .61 with vertices v$ = (0, 0), V\ = (1, 0), w 2 = (0, 1). For 
^ corresponding to the standard inner product, recall that the constant terms of the y, of the vertex cones 
are 1/2, 3/8, and 3/8. According to the discussion above, we have that the constant terms of the is's are 
found by subtracting from 1/2: 

u(K )(0) = \ K#i)(0) = \ K# 2 )(0) = \. 

One can then check that the volume of this triangle is given by summing v(F)(0) values times the number 
lattice points in the relative interior of F. Indeed, since the only lattice points are the vertices, one has 

Vol(P) = - + - + - = -. 
V ; 4 8 8 2 



SUM-INTEGRAL INTERPOLATORS AND THE EULER-MACLAURIN FORMULA FOR POLYTOPES 



21 



For the triangle from Example l5.8[ with vertices = (0, 0), V\ — (2, 0), v 2 — (0, 1), one similarly computes 
the v(K) constant terms (again using the standard inner product): 

f(*o)(0) = \ K^i)(0) = ^ v{K 2 ){Q) = ~. 

This time there is one lattice point in the relative interior of an edge. The corresponding v constant term is 
1/2, according to Proposition 16.31 Thus, we check 
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